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SECTION  I 


INTRODUCTION 


A significant  percentage  of  the  life  cycle  cost  of  a 
fleet  of  aircraft  is  due  to  the  cost  required  to  maintain  the 
aircraft  for  the  safe  performance  of  its  required  mission.  Let 
structural  cost  of  ownership  denote  that  portion  of  the  total 
life  cycle  cost  which  is  attributable  to  the  inspections  and 
repair  of  the  primary  aircraft  structure.  This  structural  cost 
of  ownership  depends  on  the  loads  environment  that  will  be 
encountered  by  Lhe  individual  aircraft  of  the  fleet  as  well  as 
damage  tolerance  properties  that  resulted  from  the  design  and 
manufacturing  processes. 

The  routine  use  of  an  operational  aircraft  results  in  the 
structure  being  subjected  to  a large  number  of  significant  load 
cycles  which  occur  randomly  in  time  and  which  have  randomly 
varying  magnitudes.  These  repeating  load  cycles  can  eventually 
weaken  the  structure  to  the  extent  that  an  applied  load  which 
is  within  the  original  design  capability  will  cause  failure  of 
the  structure.  The  fatigue  damage  resulting  from  the  applied  load 
cycles  is  currently  being  characterized  in  terms  of  the  growth 
of  cracks  at  the  critical  locations  in  the  structure.  The  crack 
growth  depends  on  the  material  and  geometry  of  the  structural 
element,  the  crack  length  at  the  time  of  a stress  cycle,  and  the 
magnitude  and  order  of  the  stress  cycles  imposed  by  the  random 
loads  on  the  aircraft.  Thus,  given  a structural  detail  with  a 
crack  of  specified  initial  length  and  the  stress  history  which 
will  be  imposed  on  the  detail,  crack  length  as  a function  of 
flight  time  can  be  calculated.  This  function  plays  a key  role 
in  the  decision  making  processes  employed  to  insure  the  structural 
integrity  of  aircraft. 
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Structural  reliability  analyses  have  been  derived  which 
present  methods  of  calculating  the  probability  the  structure  will 
not  fail  for  T flight  hours  in  an  operational  environment. 
(References  1 through  5) . These  analyses  take  into  account  some 
or  all  of  the  following  factors:  time  to  crack  initiation;  initial 
flaw  size  distribution;  random  load  sequences  representative  of 
the  real  environment;  random  effect  of  inspections;  repair  quality; 
and  material  variability.  However,  the  practical  application  of 
these  methods  require  extensive  data  bases  which,  in  general, 
are  not  available.  Thus,  although  analytical  models  are  available 
for  calculating  structural  reliability  under  different  inspections, 
they  cannot  be  applied. 

To  assure  flight  safety  of  structure,  the  Air  Force  has 
adopted  an  inspection  interval  which  is  based  on  predicted  crack 
growth  of  the  largest  flaw  in  a critical  location  which  could 
pass  the  quality  control  inspection.  To  allow  for  variability  in 
crack  growth  and  inspection  reliability,  the  structure  must  be 
inspected  within  half  the  time  interval  required  for  the  maximum 
flaw  to  grow  to  an  unstable  length.  To  date,  this  flight  safety 
criteria  has  been  the  determining  factor  in  scheduling  structural 
inspections . 

Now  during  the  periodic  inspections,  all  cracks  which  are 
detected  are  repaired  and,  in  general,  the  cost  of  repair  depends 
on  the  size  of  the  crack.  Thus,  the  structural  cost  of  ownership 
of  a fleet  of  aircraft  is  the  cost  of  all  inspections  and  repairs 
plus  the  cost  of  any  aircraft  which  may  be  lost  due  to  catastrophic 
failure.  Since  the  number  of  cracks  that  are  detected  and  repaired 
at  an  inspection  is  a random  variable  (depending  on  both  the 
distribution  of  crack  sizes  at  the  time  of  the  inspection  and  the 
probability  of  detecting  the  crack  for  the  inspection  method 
being  employed) , the  structural  cost  of  ownership  is  also  a random 
variable  and  should  be  described  in  statistical  terms. 

This  report  presents  the  results  of  a study  to  estimate 
the  structural  cost  of  ownership  due  to  the  detection  and  repair 
of  cracks.  A statistical  model  formulated  by  Yang  (Reference  6) 
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constituted  the  essence  of  the  approach.  This  model  accounted 
for  an  initial  flaw  size  distribution,  flaw  growth  by  a deter- 
ministic equation,  probability  of  flaw  detection  as  a function  of 
flaw  size,  renewal  of  repaired  fluws,  and  costs  of  inspection 
and  repair.  The  basic  model  was  generalized  to  permit  variability 
in  flaw  growth  due  to  flight  to  flight  variation  of  usage 
(Reference  7) . Since  the  model  simulates  the  growth  and  repair 
of  a distribution  of  cracks  in  a structural  detail,  the  output 
of  the  model  is  expressed  in  statistical  terms.  The  results 
cannot  be  used  to  infer  crack  size  or  cost  of  a repair  for  a 
specific  airplane  tail  number. 


Section  II  presents  a summary  of  the  method  of  calculating 
the  crack  size  distribution  at  each  inspection  interval  and  the 
expected  cost  of  inspection  and  repair.  Also  included  is  a 
discussion  on  the  particular  method  selected  for  modeling  crack 
growth  and  for  accounting  for  the  flight  to  flight  variation  of 
usage.  The  program  for  calculating  predicted  costs  is  presented 
in  Section  m with  a description  of  methods  for  estimating  the 
statistical  parameters  of  the  distributions  of  input  parameters. 
Examples  of  the  use  of  the  program  are  presented  in  Section  IV. 
The  examples  are  actually  a series  of  sensitivity  studies  to 
investigate  the  effect  of  variations  in  input  parameters  on  the 
expected  cost  of  ownership. 


SECTION  II 
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CRACK  REPAIR  COST  MODEL 

The  objective  of  the  structural  cost  of  ownership  model 
is  to  estimate  the  fleet  costs  associated  with  the  inspection 
and  repair  of  potential  cracks  at  a critical  location  on  the 
airframe.  Since  total  maintenance  costs  can  be  obtained  by  summing 
the  estimated  costs  for  the  individual  locations,  this  study  has 
been  restricted  to  a single  location.  The  following  paragraphs 
present  a summary  of  the  estimation  method.  The  basic  approach  is 
due  to  Yang  (Reference  6).  However,  the  method  of  calculating 
crack  growth  and  the  method  of  accounting  for  the  effect  of  flight 
to  flight  variability  are  different  from  Yang's  approach  and  a 
discussion  of  these  features  is  also  presented. 

2.1  MATHEMATICAL  FORMULATION 

Assume  that  a fleet  of'  aircraft  contains  N structural 
details  which  will  be  subjected  to  statistically  equivalent 
stress  environments.  There  may  be  one  structural  detail  (critical 
location)  on  each  aircraft  or  there  may  be  multiple  details  on 
each  aircraft.  The  important  elements  are  that  the  properties 
which  determine  crack  growth  are  equivalent  for  all  of  the  details. 
Since  no  two  aircraft  will  be  flown  the  same,  the  actual  stress 
histories  endured  by  the  aircraft  will  be  different.  However,  it 
is  assumed  that  stress  environments  are  statistically  equivalent 
in  their  effect  on  crack  growth. 

As  the  individual  aircraft  perform  their  operational  missions, 
cracks  will  initiate  and  grow  at  the  critical  point.  Inspections 
will  be  performed  at  increments  of  T flights  (or  an  equivalent 
number  of  flight  hours)  and,  at  each  inspection,  all  detected 
cracks  will  be  repaired.  The  cost  of  repair  will  be  a function  of 
the  size  of  the  crack  and  it  is  assumed  that  there  are  4 ranges 
of  crack  size  with  a constant  repair  cost  in  each  range.  It  is 
assumed  that  the  initial  quality  of  the  details  can  be  characterized 
by  an  equivalent  initial  flaw  size  distribution  which  depends  on 
the  initial  manufacturing  and  quality  control  processes. 
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(References  8 and  9) . Similarly,  the  repaired  details  are 
characterized  by  an  equivalent  repaired  flaw  size  distribution. 
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At  the  ith  inspection  cycle  (after  i T flights) , the  cost 
associated  with  the  inspection  and  repair  of  the  detail  under 
consideration  is  given  by 


4 

C.  = I.  + Z n (i,  j)  C (i,  j)  (1) 

j-1 

where 


1^  = cost  of  ith  inspection 

n(i,j)  = number  of  level  j cracks  discovered  and 
repaired  at  the  ith  inspection 

c(i,j)  = cost  of  repairing  one  level  j crack  at 
ith  inspection 

Since  n(i,j)  depends  on  the  distribution  of  crack  sizes  immediately 
prior  to  the  ith  inspection  and  the  probability  of  crack  detection 
for  the  inspection  method,  the  cost  of  the  maintenance  action  is 
a random  variable  and  must  be  described  in  statistical  terms. 

To  model  the  random  variable  n(i,j),  assume  that  crack 
growth  and  the  probability  of  crack  detection  are  statistically 
independent  across  structural  details.  Let  p ( i , j ) denote  the 
probability  of  detecting  (and  repairing)  a level  j crack  at  the 
ith  inspection.  Then,  n(i,j)  will  have  a binomial  distribution 
with  parameters  (N,  p(i,j)).  The  expected  value  and  variance 
of  the  number  of  detected  j level  cracks  at  inspection  i are  given 
by 


n (i , j ) = N p (i , j ) 

Sn  (i , j ) = N P (i'j)  U“P 


(2) 


Thus,  the  expected  value  and  variance  of  the  maintenance  costs 
at  the  ith  inspection  are  given  by 


F 


: 


5 


(3) 


E (ci>  = xi  + C (i,j) 


N • p(i,j) 


S2(C.)  = I C2(i,j)  N p(i,j)  [1-p  (i , j ) ] 

j = l 

and  the  expected  value  and  variance  of  the  total  maintenance  costs 
through  the  ith  inspection  (for  all  i T flights)  is  given  by 

i 

E 


E (CTOT  (i)>  = E «V 


(4) 


Var  (CTOT(i)) 


l 

E 

k=l 


(Ck> 


If  N p(i,j)  and  N (l-p(i,j))  are  both  greater  than  5,  then  the 
normal  distribution  may  be  used  to  approximate  the  binomial  and 
the  distribution  of  costs  at  the  ith  inspection  and  for  the  total 
of  all  inspections  can  be  characterized  by  the  normal  distribution 
with  parameters  as  specified  in  (3)  or  (4) . 

The  calculation  of  p(i,j)  depends  on  the  distribution  of 
crack  sizes  at  the  ith  inspection  and  the  probability  of  detecting 
a crack  as  a function  of  crack  size  for  the  particular  inspection 
technique  being  used.  Let  f^  (a)  represent  the  density  function 
of  crack  sizes  immediately  before  the  ith  inspection.  Let  H(a) 
denote  the  probability  of  detecting  a crack  of  length  a.  Then 

a . 


P (i,j)  = f 


H(y)  fi  (y)  dy 


(5) 


where  the  a ^ , j = 0,1, 2, 3, 4 define  the  crack  size  levels  and  aQ=0, 
a4=  - . 

Immediately  after  the  repair  of  the  cracks  the  density 
function  of  the  crack  sizes  is  a mixture  of  the  density  function 
of  the  cracks  which  were  not  found  and  the  equivalent  repaired 
crack  size  density  function.  If  ft  (a)  represents  the  density 
function  of  crack  sizes  immediately  after  repair,  then 


f.  (a)  = P (i)  f (a)  + tl-H(a)]  f.  (a) 
1 R 1 


(6) 


MtllU U 


where 


and 


f R <a) 


P(i)  = 


= equivalent  repaired  crack  size  density 
function 


proportion  of  all  cracks  which  were  repaired 
4 

I P (i,j) 

3=1 

(7) 


fi  (y)  dy 


Using  equations  (4)  through  (7) , the  expected  value  and 
variance  of  the  inspection  and  repair  costs  can  be  obtained  at 
each  inspection  period  and  the  distribution  of  the  crack  sizes 
immediately  after  the  inspection  can  be  calculated.  The  only 
remaining  element  which  must  be  modeled  is  the  density  function 
of  the  crack  sizes  immediately  before  the  inspection.  This 
distribution  is  obtained  iteratively  bv  assuming  cracks  grow 
in  accordance  with  the  equation 


a (T) 


b x (T) 
e 


(8) 


where  a^  = crack  length  after  ith  inspection  and  repair 

b = constant  depending  on  material  and  detail  geometry 

x (T)  = normally  distributed  random  variable  with  parameters 
and  oT  which  accounts  for  the  random  stress 
environment  encountered  by  a particular  detail 
during  the  T flights. 


This  crack  growth  model  is  discussed  in  detail  in  paragraph  2.2. 

Given  the  distribution  of  flaw  sizes  immediately  following 
an  inspection  and  repair,  fT  (a) , and  the  crack  growth  equation, 
(8) , the  density  function  of  flaw  sizes  immediately  before  the 
next  inspection  (after  T flights)  can  be  obtained  as  follows. 
First  the  cumulative  distribution  of  flaw  sizes  aftt  .*■  T flights 
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is  calculated 


Fi+i(a) 


= P U1+i  < a> 


= P {ai  exp  (b  x ) < a} 


- a e 

■f  f 

_ m — O 


a exp  (-b  x) 


^(yJdy  g(x)  ax  (9) 


where  g (x)  is  the  normal  density  function 


g (x)  = 


/2rr  ar 


exp  - [(x-ut)2/2 a2] 


The  density  function  is  then  obtained  as 


fi+l  (a)  = 


■Si-  Fi+1  (a) 


-I 


exp(-b  x)  f ^ (a  exp(-b  x)  ) g(x)dx  (11) 


Further,  letting 

z = x-yT/oT 

then,  for  all  practical  purposes 

/*  3 


f i+1 (a) 


■/ 


exp [-b  (yT  + 20t)]  (a  exp  [-b  (yT+ zjt) ] ) g (z) dz 


with 


g(z)  = 


/2tt 


- (Z  2/2) 


The  above  mathematical  model  for  simulating  crack  growth, 
inspection,  and  cost  of  repair  has  been  programmed  to  provide 
output  at  each  of  the  pre-specif ied  inspection  cycles  through  a 
pre-defined  lifetime.  The  simulation  process  starts  with  an 
equivalent  initial  crack  size  distribution  from  which  the  distribution 


of  flaw  sizes  immediately  before  the  first  inspection  is  calculated, 
equation  (12)  . The  costs  of  inspection  and  repair  for  this  first 
inspection  cycle  are  calculated  using  equations  (5)  and  (3) . 

(At  subsequent  inspections,  total  costs  are  also  calculated 
using  equation  (4).)  Equations  (7)  and  (6)  then  provide  the 
density  function  of  flaw  sizes  immediately  after  the  first 
inspection  and  this  density  is  used  to  initiate  the  second  cycle 
of  simulated  operation.  The  process  is  repeated  for  a specified 
number  of  cycles  (aircraft  life) . This  process  is  described  in 
greater  detail  in  Section  III,  which  presents  a description  of  the 
computer  program  and  its  use. 

2.2  DISCUSSION 

The  mathematical  model  presented  in  the  previous  paragraph 
depends  on  basic  assumptions  regarding  the  process  of  crack  growth 
and  repair  during  operational  usage.  The  assumption  regarding 
independence  of  crack  detection  between  details  is  a probabilistic 
assumption  which  is  non-controversial  and  permits  characterization 
of  the  cost  of  repairs  in  terms  of  the  binomial  distribution.  A 
second  class  of  assumptions  are  those  which  relate  to  the  equivalent 
distribution  of  crack  sizes.  This  method  of  characterizing  initial 
and  repair  quality  has  met  with  general  acceptance  in  the  Air 
Force  and  the  only  difficulty  in  the  application  is  the  specifica- 
tion of  the  proper  distributional  family  and  the  values  of  the 
parameters  of  the  distribution.  Similarly,  the  concept  of  the 
probability  of  detection  of  a crack  as  a function  of  crack  size 
is  well  accepted  but  a specific  distribution  form  and  specific 
parameter  values  for  an  NDI  method  may  not  be  known.  The  mathe- 
matics of  the  crack  repair  cost  model,  however,  do  not  depend  on 
the  specific  family  or  its  parameters  for  either  of  the  latter 
types  of  distributions. 

The  method  of  simulating  crack  growth  in  equation  (8)  is. 


however,  a different  concept  which  requires  discussion.  In  the 
following  two  paragraphs,  a heuristic  argument  is  presented  for 
the  use  of  this  model. 


2.2.1  Crack  Growth 


The  crack  growth  model  of  Yang  is  based  on  the 
assumption  that  incremental  crack  growth  during  the  ith  flight 
is  given  by 

Aai  = b a?..!  r c * 1 (13) 

After  t flights,  a crack  of  initial  length  aQ  would  be  given  by 

a = aQ  [1  + (1-c)  b t aoC_1]  1/1'c  (14) 

In  particular,  Yang  states  that  for  aluminums,  c = 2,  so  that 


1-b  t a 

o 

This  model  is  completely  deterministic  in  that  a given  initial 
crack  length  completely  determines  the  history  of  the  growth  of 
the  crack. 

Now  crack  growth  tests  from  simulated  operational 
usage  have  indicated  that  a model  of  the  form 

a = aQ  exp  (b  t)  (16) 

also  fits  the  data  reasonably  well  for  values  of  t less  than  the 
design  life.  As  an  example.  Figure  1 presents  crack  length  data 
(with  a least  squares  fit  of  the  exponential  model  superimposed) 
from  crack  growth  tests  during  the  F-5  damage  tolerance  assessment. 
Reference  (10) . This  example  is  typical  of  the  early  stages  of 
crack  growth  where  the  growth  is  very  slow.  Thus,  there  is  evidence 
to  indicate  that  the  exponential  model  may  provide  a reasonable 
(if  not  good)  approximation  for  crack  length  as  a function  of  time. 
This  is  particularly  true  for  the  earlier  stages  of  crack  growth 
(up  to  half  crack  life)  which  are  the  periods  of  major  interest  in 
the  calculation  of  the  expected  cost  of  crack  repair. 

The  exponential  model  results  from  the  assumption 
that  the  crack  growth  during  flight  i is  directly  proportional  to 
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rack  Lengt 


the  crack  length  at  the  start  of  the  flight: 

Aa^  = b ai-:l  (17) 

Thus,  the  exponential  model  is  the  special  case  of  equation  (13) 
with  c = 1.  Given  a crack -length  versus  number- flights  data  set 
for  a detail  of  interest  under  the  anticipated  loading  history, 
the  parameter  of  the  model  can  easily  be  calculated,  e.g.,by  a 
least  squares  estimate. 

To  consider  the  crack  growth  model  from  another 
approach,  assume  that 

= b (AK)m  (18) 

where 

AK  = 8 (a)  / Tra  Ao  (19) 

8 (a)  = factor  dependent  on  crack  length  and  geometry. 

This  crack  geometry  model  has  often  fit  empirical  data  in  the  mid 
ranges  of  crack  growth  and  provides  an  acceptable  fracture  mechanics 
model  for  fleet  modeling  where  actual  stress  time  histories  in  each 
aircraft  will  not  be  available  (Reference  11) . For  aluminum, 
values  of  m in  the  range  of  3 to  4 are  commonly  observed. 

During  flight  i an  aircraft  will  experience  n^  load 
cycles  and  the  crack  will  grow  an  increment  Aa^  where 


Aa.  = 


A a . . 


Assuming  negligible  error  results  if  crack  length  at  start  of 
flight  is  used  to  calculate  all  increments  during  the  flight, 
substituting  from  (18)  and  (19)  fields 


Aa.  = E b [/F  8(a.)/aT  Ao.]1 

j = l 1 3 


= b \\m/i  [8m  (a.)  a.  m/z]  E1  (Ao.)m 

11  j=l  3 


Since  the  exponential  model  has  been  observed  to  fit  crack  length 
data  in  some  applications,  it  may  be  reasonable  to  assume  that  for 
the  Aa  ranges  of  interest  that 

fy.  m t * in/  2 

B (a^)  ai  ' 2 a. 

Thus,  there  exists  a constant  b such  that  the  crack  growth  during 
the  ith  flight  is  given  by 


where 


Aa.  = b a X. 
i l 


X.  = I (Ac.)' 
j = l D 


Since  the  number  and  magnitude  of  the  stress  cycles  to  be  encountered 
in  a randomly  selected  flight  are  not  known,  Xi  can  be  considered 
as  a random  variable  which  accounts  for  the  stress  cycles  of  the 
flight.  If  the  crack  growth  tests  for  the  structural  detail  were 
based  on  a random  spectrum  which  meets  some  baseline  usage,  the 
exponential  crack  growth  curves  could  be  considered  as  resulting 
from  the  statistical  model 


E (Aai)  = b E(Xi)  a.  (23) 

The  statistical  approach  introduced  by  equation  (23)  will  be  further 
developed  in  the  next  paragraph. 

2.2.2  Flight  by  Flight  Variability 


In  the  previous  paragraph,  a heuristic  argument  was 
presented  for  modeling  the  incremental  crack  growth  during  flight  i 
by  the  equation 

t\aL  = b ai_1  Xi  (24) 

where  b is  a constant  determined  from  specimen  test,  ^ is  the 
crack  length  after  i-1  flights  and  x^,  defined  by  equation  (23), 
is  a random  variable  which  characterizes  the  non-deterministic 
stress  history  to  be  encountered  during  the  ith  flight.  Using 
this  relation,  it  is  easily  shown  that  the  crack  length  after  t 
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flights  is  given  by 


from  whence 


(1+b  Xi) 


£n  a.  - in  a 
t o 


t 

Z £n  (1+b  X. ) 
i=l  1 


(25) 


Since  b X^  = Aa^/a^_^  <<|  over  the  crack  range  of  interest,  it 

can  be  shown  that 


or 


where 


Kn  afc 


£n  a 

o 


t 

b Z X. 
i=l  1 


t 

at  = ao  exp  (b  E xi> 
i=l 


bX  (t) 


t 

X(t)  = Z X.. 

i=l  1 


(26) 


By  the  Central  Limit  Theorem,  when  t is  large  (>50)  , X(t)  can  be 
assumed  to  have  a normal  distribution  with  mean  and  variance  given  by 


“t  * 1 E(V 

2 

ofc  = t var  (X^ 


(27) 


Thus,  it  can  be  assumed  that  the  random  usage  factor  of  the  crack 
growth  model  has  a normal  distribution  as  specified  in  equation  (10) . 

Now  the  crack  growth  rate  parameter,  b,  will  be 
estimated  from  a baseline  spectrum  which  represents  the  anticipated 
stress  environment.  Rather  than  calculate  E(X.)  for  this  spectrum 
to  obtain  a growth  rate  parameter  that  is  independent  of  the  baseline 
spectrum,  a single  growth  rate  parameter  which  includes  the  baseline 
environment  can  be  calculated  and  the  flight  by  flight  variation  can 
be  expressed  in  terms  of  the  coefficient  of  variation.  Thus,  if  b is 
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calculated  from  the  baseline  spectrum  and  the  aircraft  of  the  fleet 
are  flying  this  spectrum,  then 


I 


°v 

a - /t  • i 

E(Xi) 

In  the  absence  of  specific  stress  history  data,  neither  E(X^)  or 
ov.  can  be  known.  However,  coefficients  of  variation  of  the  Miner's 
damage/ flight  have  been  calculated  in  a few  cases  and  these  values 
provide  an  indication  of  the  relative  magnitude  of  the  variability 
of  crack  growth  per  flight.  These  values  range  from  about  unity 
(o=u)  for  specific  well  defined  missions  to  approximately  4 (o=4u) 
for  composites  of  several  missions  of  a fighter/attack  aircraft. 

Such  large  coefficients  of  variations  would  also  be  expected  for  the 
random  variable  X^,  equation  (23),  since  both  metrics  are  a function 
of  the  number  and  magnitude  of  the  same  random  number  of  random 
stress  cycles. 

Since  the  crack  growth  model  is  directly  dependent 
upon  E(X^),  equation  (27)  suggests  a simple  test  of  equality  of 
sample  means  as  a method  for  determining  if  operational  usage  is 
different  from  that  of  the  baseline  spectrum.  Sample  size  considera- 
tions could  also  be  expressed  in  terms  of  number  of  flights  to  be 
monitored  to  obtain  a pre-determined  precision  in  the  value  of 
E(X^).  These  statistical  considerations  have  been  considered  in 
Reference  (12)  which  used  average  Miners  damage/flight  as  the 
monitoring  parameter.  However,  the  statistical  considerations  would 
transfer  directly  to  the  new  parameter  E(X^). 

Note  that  the  model  for  accounting  for  flight  to 
flight  variability  in  crack  growth  during  operational  usage  depends 
on  the  assumption  that  the  stress  environment  encountered  by  a 
particular  aircraft  during  a flight  is  a random  selection  of  all 
such  flights  of  all  aircraft  of  the  fleet.  If  a particular  aircraft 
flies  a snectrum  that  is  biased  toward  the  severe  side,  it  will,  of 
course,  experience  crack  growth  greater  than  that  predicted  by  the 
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model.  Individual  aircraft  tracking  is  required  to  account  for 
specific  aircraft  usage. 
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SECTION  III 


APPLICATION  OF  MODEL 


A computer  program  was  written  to  calculate  the  predicted 
crack  repair  costs  (PCRC)  in  accordance  with  the  model  expressed 
in  the  previous  paragraphs.  Figure  2 presents  a schematic  of  the 
calculations  and  indicates  the  iterative  nature  of  the  computation. 
The  issue  specifies  a set  of  input  data  related  to  flaw  sizes, 
cost  of  inspections  and  repairs,  and  number  of  inspection  cycles. 

At  each  inspection,  the  program  outputs  the  percentage  of  cracks 
in  each  repair  interval  which  was  identified  by  the  inspection 
process,  the  total  percentage  of  cracks  found,  the  expected  value 
and  standard  deviation  of  the  inspection  and  repair  cycles,  and 
the  cumulative  cost  totals.  At  the  users  option,  the  program  will 
also  print  out  the  density  functions  of  the  crack  sizes  immediately 
before  and  after  the  inspections. 

In  writing  the  program,  specific  assumptions  regarding  the 
equational  form  of  distributions  had  to  be  made.  In  the  following 
paragraphs,  all  of  the  input  parameters  required  by  the  program 
are  discussed.  The  input  parameters  are  considered  in  three 
groupings:  (1)  flaw  size  parameters,  (2)  inspection  and  repair 
parameters,  and  (3)  program  parameters.  The  last  paragraph  of 
this  section  presents  the  program  and  its  user's  manual. 

3.1  FLAW  SIZE  PARAMETERS 

The  input  parameters  related  to  flaw  size  are  the  equivalent 
initial  and  repaired  flaw  size  distributions  and  the  crack  growth 
parameters.  Each  of  these  is  discussed  in  the  following. 

3.1.1  Initial  Flaw  Size  Distributions 

The  equivalent  initial  flaw  size  distribution  is 
a concept  used  to  quantify  the  concept  of  manufacturing  quality. 
Direct  measurements  of  initial  flaw  sizes  cannot  be  made  so  that 
the  distributional  family  that  is  used  to  represent  initial  quality 
must  be  inferred  from  other  data.  One  method  is  to  measure  the 
cracks  in  a structural  item  after  exposure  to  a stress  environment 
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and  using  that  stress  environment  to  back  calculate  the  "initial" 
crack  size  at  time  zero.  Reference  8 presents  the  distributions 
obtained  from  A-7D  and  F-4  C/D  aircraft.  The  A-7D  initial  hole 
quality  data  is  represented  reasonably  well  by  the  Weibull  family 
of  distributions.  The  F-4  C/D  data  deviate  from  the  Weibull  family 
particularly  for  the  very  small  crack  sizes.  (The  F-4  C/D  have 
been  fit  by  a Johnson  Su  model  but  the  use  of  this  particular  4 
parameter  family  cannot  be  justified  for  any  other  data  set. 

Further,  there  is  some  question  regarding  the  homogeneity  of  the 
data  points  of  the  F-4  C/D  data  set  (Reference  7).) 

Since  the  two  parameter  Weibull  family  of  distri- 
butions provides  four  different  slopes  of  the  density  function,  since 
it  is  a commonly  accepted  and  used  distribution,  since  it  is 
analytically  easy  to  manipulate,  and  since  it  is  a reasonable 
model  for  at  least  one  data  set  (the  A-7D) , this  family  was 
programmed  in  the  model  for  the  specification  of  equivalent 
initial  hole  quality.  If  the  user  prefers  another  two-parameter 
distribution  a new  FORTRAN  equation  must  be  entered  for  the  FUNCTION 
FO (X)  subroutine  with  the  parameters  defined  as  FOALPH  and  FOBETA. 

If  a distribution  with  three  or  more  parameters  is  desired, 
additional  changes  to  the  input  logic  will  be  required. 

The  density  function  and  cumulative  distribution 
function  of  the  two-parameter  Weibull  distribution  are  given  by 

f (y)  = | [y/8]01""1  exp  - [y/8]a  (29) 

F (y)  = 1 - exp  - [y/B]a 

The  shape  of  the  distribution  is  determined  from  a while  3 is  a 
scale  parameter  and  63  percent  of  the  distribution  is  always  less 
than  8 (i.e.,F(8)  = 0.63).  8 is  also  known  as  the  characteristic 

value.  Figure  3 presents  plots  of  the  Weibull  density  functions 
for  5 values  of  the  shape  parameter.  For  each  value  of  a,  the 
scale  parameter,  8/  was  determined  such  that 

F (0.050)  = 0.999 
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Weibuil  Density  Functions  With  99.9  percent  Less  Than  0.050  in. 


This  vlaue  was  selected  since  it  is  representative  of  the  type  of 
equivalent  initial  quality  that  is  required  by  Mil  Std  1530A. 

Note  that  for  ct<l,  the  density  function  is  asymptotic  to  the 
vertical  axis  and  a high  percentage  of  the  cracks  is  very  small. 

For  a=l,  the  exponential  distribution  results.  For  a>l,  the  density 
function  will  have  1 or  2 points  of  inflection  and  as  a increases, 
the  scatter  about  the  median  decreases  and  the  median  increases. 

The  data  from  the  A-7D  and  F-4  C/D  initial  hole  quality  studies 
result  in  values  of  a between  2 and  3. 

Given  a set  of  equivalent  initial  hole  quality 
data  which  is  applicable  to  the  structural  detail  of  interest, 
the  Weibull  parameters  can  be  estimated  through  graphical  (Weibull 
probability  paper)  or  analytical  (maximum  likelihood)  methods. 

Given  the  data  points  the  appropriateness  of  the  Weibull  family 
could  also  be  tested.  In  the  absence  of  data,  however,  the  values 
of  the  parameters  must  be  inferred  from  other  information.  For 
example,  if  a quality  control  process  has  been  shown  to  be  effective 
at  detecting  flaws  of  a given  magnitude  with  a high  probability, 
this  information  can  be  interpreted  as  a point  on  the  cumulative 
function. 

In  general,  either  a or  3 can  be  calculated  from 
a specified  value  of  the  other  and  from  one  specified  point  on 
the  cumulative  distribution  functions.  Given  a and  the  point 
F(Y  ) = P,  then 

*r 

3 = ^2 (30) 

(-£n  (1-p)  ] 1/cx 

For  example,  if  it  is  assumed  that  the  shape  parameter,  a,  is 
two  in  accordance  with  the  A-7D  data  and  it  is  shown  that  99.9 
percent  of  all  initial  flaws  are  less  than  50  mills,  i.e.,  flaws 
greater  than  50  mills  are  eliminated  by  quality  control,  then 

F ( 0 . 050)  = 0.999 
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and 


6 


0.050 

[ -dn(l-0.999)  ] 1/2 


= 0.019  inches 

Given  g and  the  point  F(yp)  = P,  then 


_ dn  [-dn  (1-p) ] 
dn  yp  - in  B 


(31) 


For  example,  if  assumed  as  before  that  99.9  percent  of  all  initial 
flaws  are  less  than  50  mills  and  that  the  characteristic  value  of 
an  initial  flaw  size  is  20  mills,  then 


a=  dn [-dn  (1-0.999)] 
dn (0. 050)- dn (0.020) 

= 2.11 

If  two  points  on  the  cumulative  distribution 

function  are  specified  (y  , P ) and  (y  , P ) then 

11  2 2 


where 


dn  y -dn  Y 

B=exp  ( l- 7- — ] 

A-l 


(32) 


dn  [-dn(l-P) ] 

A = 2- — 

dn [-dn (1-P  ) ) 
i 

and 


dn [-dn (1-P  ) ] 

a= 1 (33) 

dn  Y - dn  6 
1 

For  example,  if  99.9  percent  of  initial  flaws  are  less  than  50  mills 
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and  50  percent  of  the  initial  flaws  are  less  than  20  mills,  then 


A = &n  (0.001)]  _ _5  # 27  306 

fin  [-Jin  (0. 50)  ] 

g A In  (0. 020) -fin  (0.050)  ] 
A- 1 

= 0.0231 


a = fin  [-fin  (1-0.50)] 

in  ( 0 . 020) -in  (0.0231) 

= 2.51 

In  the  program,  the  shape  and  scale  parameters 
of  the  initial  crack  quality  distribution  are  denoted  by  FOALPH 
and  FOBETA,  respectively. 

3.1.2  Crack  Growth  Parameters 

There  are  four  parameters  in  the  program  which 
are  identified  as  crack  growth  parameters.  These  are  (1)  the 
number  of  flights  between  inspections,  UT;  (2)  the  average  crack 
growth  rate  constant  for  the  baseline  usage,  UB;  (3)  the  bias 
factor  which  reflects  change  in  usage  from  baseline,  UM;  and 
(4)  the  flight  by  flight  variability  factor  which  is  the  coefficient 
of  variation  of  the  flight  by  flight  measure  of  the  stress  environ- 
ment, US. 

The  number  of  flights  between  inspections  is 
entirely  at  the  discretion  of  the  user.  If  inspection  intervals 
are  planned  in  terms  of  flight  hours  rather  than  number  of  flights, 
the  average  flight  length  provides  an  acceptable  conversion  to 
number  of  flights.  This  parameter  is  considered  as  a crack  growth 
parameter  since  it  is  required  in  the  statistical  distribution 
which  describes  effect  of  flight  to  flight  variability  on  crack 
growth,  equations  (28)  and  (12) . 
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The  crack  growth  rate  parameter  b,  or  UB  in  the 
program,  can  be  estimated  from  the  results  of  random,  flight  by 
flight  crack  growth  tests.  A least  squares  fit  of  equation  (16) 
would  provide  an  estimate  of  b and  if  multiple  specimens  were 
tested  to  the  same  spectrum,  an  average  of  b values  would  be 
used  in  the  model.  Similarly,  a crack  length  vs.  number  of  flights 
curve  could  be  generated  by  a cycle  by  cycle  crack  growth  equation 
for  the  baseline  spectrum  and  b could  be  determined  to  provide  a 
fit  of  equation  (16)  to  this  curve. 

In  the  absence  of  data  or  for  analytical  studies 
not  related  to  specific  details,  reasonable  values  of  b can  be 
inferred  from  assumptions  regarding  the  expected  safe  life  of 
a detail  and  the  initial  quality  distribution  of  a detail.  For 
example,  if  it  has  been  demonstrated  or  assumed  that  flaws  of  a 
magnitude  greater  than  aQ  will  be  detected  by  the  manufacturer, 
that  ac  is  the  critical  crack  length  for  the  detail,  and  that  the 
life  of  a detail  is  Tc  flights,  then 


and 


a = a e 
c o 


bTc 


b 


Jin 


iln  ap 


(34) 


If  the  safe  life  estimate  T is  conservative  (short)  then  the  crack 

c 

growth  parameter  will  be  conservative  (large) . 


To  provide  an  indication  of  the  magnitude  of  the 
crack  growth  rate  parameter , Table  I presents  the  least  squares 
estimates  obtained  from  specimen  crack  growth  test  for  the  F-5  E/F 
ai-rcraft.  Reference  10.  Two  specimens  were  tested  under  identical 
test  conditions  so  that  the  observed  differences  in  b values  for 
the  duplicates  are  due  to  uncontrolled  sources  of  variation. 

Note  that  each  individual  crack  growth  curve  was  fit  very  well  by 
the  exponential  model  but  that  for  some  specimen  types,  signifi- 
cant differences  in  estimated  b values  were  obtained  from  duplicate 
tests . 


24 


The  usage  bias  and  flight  by  flight  variability 
parameters  can  be  considered  together.  If  the  anticipated  stress 
environment  is  equivalent  to  that  of  the  baseline  spectrum,  there 
will  be  no  usage  bias  and  the  usage  bias  factor  is  unity.  If  the 
anticipated  environment  is  different  from  the  baseline,  the  ratio 
of  anticipated  usage  severity  to  baseline  is  the  usage  bias 
parameter.  Equation  (23)  expressing  a flight  by  flight  usage 
parameter  in  terms  of  the  summation  of  a power  of  stress  changes 
is  a parameter  which  may  provide  a measure  of  this  usage.  Further 
development  is  required  on  the  application  of  this  parameter. 

TABLE  I 

Crack  Growth  Parameters  From 
F-5  E/F  Specimen  Tests 


Specimen 

Identification 

Max  Spectrum 
Stress 

■ r 

/\ 

ao 

b (X  104) 

01  W 1 

32910 

.020 

2.93 

01  W 2 

32910 

.018 

2.40 

02  W 1 

29120 

. 019 

2.07 

02  W 2 

29120 

.024 

1.68 

03  W 1 

30150 

.029 

2.62 

03  W 2 

30150 

.029 

2.20 

04  W 1 

33730 

. 027 

2.62 

04  W 2 

33730 

. 022 

1.56 

05  W 1 

30810 

.037 

2.77 

05  W 2 

30810 

.036 

4.38 

06  W 1 

± 29030 

.054 

1.12 

06  W 2 

± 29030 

.049 

1.19 

OS  W 1 

36880 

.023 

4.63 

08  W 2 

36880 

.017 

5.35 

The  flight  by  flight  variability  parameter  required 
by  the  program  is  the  coefficient  of  variations  of  the  flight  by 
flight  usage  parameter.  No  experience  is  available  for  the  usage 
parameter  defined  by  equation  (23)  but  Minors  damage  calculations 
indicate  that  any  realistic  metric  of  flight  by  flight  structural 
damage  will, in  general,  have  an  erratic  distribution  with  a 
standard  deviation  at  least  equal  to  the  mean.  For  Miners  damage 
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and  some  mission  types  on  fighter  aircraft,  coefficients  of  variation 
greater  than  3 have  been  observed,  c.f.  Reference  12.  Actually,  in 
the  examples  of  Section  IV  it  will  be  shown  that  the  expected  cost 
of  crack  repair  is  extremely  insensitive  to  this  parameter.  This 
point  will  be  discussed  later. 

3.1.3  Repaired  Flaw  Size  Distribution 

After  a crack  has  been  detected  and  repaired,  it 
is  assumed  that  the  detail  has  an  equivalent  crack  that  is  described 
by  the  repaired  flaw  size  distribution.  Thus,  the  exact  crack 
size  for  a particular  detail  is  not  known;  but  knowing  the  distri- 
bution would  permit  calculation  of  the  probability  that  the 
equivalent  crack  is  in  a specific  interval.  Further,  if  a large 
group  of  details  is  repaired,  the  repaired  flaw  size  distribution 
would  provide  the  percentage  of  the  details  that  have  equivalent 
cracks  in  a specified  interval. 

In  the  program,  it  is  assumed  that  the  equivalent 
repaired  flaw  size  distribution  is  from  the  two-parameter  Weibull 
family.  The  reasons  for  choosing  the  Weibull  family  are  the  same 
as  those  for  selecting  the  Weibull  family  for  the  equivalent  initial 
flaw  size  distribution.  If  a different  two-parameter  family  is 
desired,  the  change  can  be  accomplished  by  making  the  appropriate 
changes  in  the  FUNCTION  FR(X)  subroutine.  The  methods  for 
estimating  the  parameters  of  the  repaired  flaw  size  distribution 
are  the  same  as  those  presented  for  the  initial  flaw  size  distri- 
bution in  Paragraph  3.1.1.  Note  that  the  model  assumes  a constant 
repair  quality  for  each  inspection  and  for  the  repair  of  all  crack 
sizes . 

3.2  INSPECTION  AND  REPAIR  PARAMETERS 

The  input  data  which  are  identified  as  inspection  and  repair 
parameters  are  the  cost  of  each  inspection,  the  function  which 
specifies  the  probability  of  crack  detection,  and  the  costs  of 
repair  for  the  details  identified  as  having  cracks. 


3.2.1 


Cob t of  Inspection 


'trUftlMItffti 


The  model  uses  a predetermined  cost  of  inspection 
at  each  inspection  cycle.  This  cost  is  independent  of  crack  size 
and  represents  a particular  detail's  share  of  the  total  cost  of 
preparing  the  airplane  for  inspection  as  well  as  the  cost  of  the 
time  required  to  inspect  the  actual  detail.  Since  the  location 
of  the  critical  detail  could  greatly  influence  the  time  required 
for  its  inspection,  the  estimate  of  the  inspection  cost  can  only 
be  based  on  an  analysis  of  the  total  inspection  plan.  For  example, 
an  accessible  critical  location  may  be  easy  to  inspect  in  the  field 
while  an  inaccessible  critical  location  may  be  inspected  only  during 
a PDM  at  a time  when  many  other  maintenance  actions  are  performed. 


The  program  is  designed  to  permit  a different  cost 
of  inspection  at  each  cycle.  This  feature  could  be  used  to  account 
for  different  inspections  occurring  in  the  field  or  in  the  depot. 

In  order  to  set  up  an  array  for  multiple  inspection  costs,  the 
user  is  asked  to  specify  whether  or  not  multiple  inspection  costs 
will  be  specified.  If  not,  the  user  enters  the  single  value.  If 
so,  the  user  enters  one  value  for  each  inspection  cycle. 

3.2.2  Probability  of  Crack  Detection 

The  non-destructive  inspection  teqhniques  currently 
in  use  in  the  Air  Force  are  not  capable  of  detecting  all  cracks 
in  the  zone  of  inspection.  Since  the  likelihood  of  a crack  being 
detected  depends  on  a large  number  of  factors  (length,  shape, 
location  of  crack;  NDI  method;  ability  of  inspector;  location  of 
inspection;  time  of  day;  etc.)  crack  detection  capability  has 
been  characterized  in  terms  of  the  probability  of  detection  (POD) 
as  a function  of  crack  length.  The  other  factors  which  influence 
crack  detection  are  accounted  for  either  by  degrading  the  POD  for 
a particular  crack  length  due  to  the  scatter  introduced  by  the 
factors  or  by  determining  POD  curves  for  major  controllable 
stratifications  of  the  population  of  all  NDI  inspections.  Thus, 
there  may  be  a POD  curve  for  each  combination  of  inspection  techni- 
que and  structure  type  but  the  personal  factors  which  influence 
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the  inspector  (time  of  day,  day  of  week,  mental  attitude,  etc) 
would  have  to  be  measured  in  terms  of  scatter  in  POD  for  a crack 
of  fixed  length. 

For  the  computer  program  of  the  model  it  was 
assumed  that  the  probability  of  crack  detection  as  a function  of 
crack  length  can  be  described  by  a cumulative  Weibull  distribution. 
This  family  of  distributions  was  selected  since  it  includes  the 
exponential  distribution  which  has  been  used  in  the  past  as  a 
model  (Reference  6) , and  also  since  it  was  used  to  characterize 
POD  in  a recently  completed  study  for  the  Air  Force  (Reference  13)  . 
Further,  due  to  the  extreme  scatter  in  the  percentage  of  inspectors 
who  detect  cracks  of  almost  equal  length,  there  are  many  families 
of  distributions  that  could  not  be  rejected  as  models  for  the  data. 
Since  the  Weibull  family  has  an  analytically  simple  form,  it  provides 
a convenient  and  acceptable  model.  If  the  user  desires, any  other 
two-parameter  family  can  be  substituted  by  making  changes  to  the 
FUNCTION  H ( X)  subroutine.  It  should  be  noted,  however,  that  H(X) 
must  approach  one  as  X increases  to  insure  the  consistency  of  the 
mathematics  of  the  model.  Probability  of  detection  curves  are 
often  shown  as  being  asymptotic  to  a value  less  than  one.  If  such 
a curve  is  desired,  it  must  be  set  equal  to  unity,  at  some  crack 
length.  The  incorporation  of  such  a POD  would  require  significant 
changes  in  the  programmed  model. 

The  probability  of  crack  detection  in  the  program 
input  requires  the  specification  of  the  shape  (RALPH)  and  scale 
(HBETA)  parameters.  These  values  should  be  estimated  for  the  type 
of  structure  and  NDI  method  that  will  be  used  during  the  inspections. 
Figure  4 presents  the  probability  of  detection  for  four  values  of 
the  shape  parameter  and  a common  characteristic  value  of  0.1  in. 

The  data  of  Reference  13  indicate  that  a characteristic  value  of 
0.1  in.  would  be  reasonable  for  eddy  current  surface  scans  around 
countersunk  fasteners  in  skin  and  stringer  wing  segments.  This 
data  also  indicates  that  the  smaller  shape  parameters  (long  tail) 
are  more  descriptive  of  current  capabilities  than  the  less  spread 
distributions  which  result  from  larger  shape  parameters. 
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3.2.3  Cost  of  Repair 

Three  types  of  data  are  required  in  the  input 
related  to  cost  of  repair:  the  number  of  details  in  the  population, 
the  crack  sizes  which  determine  the  four  levels  of  cost  of  repair, 
and  a cost  value  for  each  of  the  repair  levels. 

The  number  of  details  in  the  population  is 
required  in  the  calculation  of  the  expected  value  and  standard 
deviation  of  the  costs  of  repair  at  an  inspection.  This  is  the 
total  number  of  details  in  the  fleet  which  are  being  subjected 
to  an  equivalent  stress  environment  with  equivalent  crack 
detection  probability  and  inspection  and  repair  costs.  In  the 
simple  case  of  one  critical  location  per  aircraft,  this  parameter 
is  the  number  of  aircraft  in  the  fleet.  In  the  program,  the 
parameter  is  entered  as  the  number  of  items,  N. 

The  calculation  of  the  expected  cost  of  repair  allows 
for  four  ranges  of  crack  size  with  a (possibly)  different  repair 
cost  for  a crack  in  each  range.  This  provision  is  based  on  the 
assumption  that  it  may  be  more  expensive  to  repair  large  cracks 
than  small  cracks.  For  example,  oversizing  a fastener  hole  may 
be  an  inexpensive  and  viable  method  for  repairing  a small  crack 
but  if  a crack  exceeds  a certain  length,  an  entire  panel  may 
require  replacement.  Further,  if  the  replacement  cost  of  the 
aircraft  is  entered  as  the  cost  of  repair  for  the  largest  critical 
flaw  size,  the  expected  costs  of  catastrophic  failure  are  also 
included  in  the  model. 

There  are  four  critical  flaw  sizes  which  define  the 
four  ranges  of  cost  repair.  In  the  input  to  the  program,  these 
are  denoted  as  AA(I),  1=1,4.  The  four  ranges  are:  AA(1)  to  AA(2); 
to  AA(3) ; AA ( 3 ) to  AA(4) ; and  greater  than  AA(4).  If  any  crack 
that  is  discovered  is  repaired,  then  AA(1)=0.  In  the  case  of  cracks 
from  fastener  holes  typical  values  for  AA(2),  AA(3),  and  AA(4)  would 
be  determined  such  that  repair  is  accomplished  by  oversize 
drilling  and  installing  a new  fastener.  For  example  AA(2)=0.016, 
AA(3)=0.032,  and  AA(4)=0.064.  For  this  set  it  would  be  assumed  that 
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replacement  (maximum  repair)  is  required  for  cracks  longer  than 
0.064. 

The  cost  to  repair  a crack  in  each  of  the  ranges 
would  depend  on  the  location  and  would  be  estimated  from  current 
experience.  It  is  possible  that  different  types  of  repair  could 
be  performed  at  different  costs  and  the  result  would  be  different 
equivalent  repaired  crack  size  distributions.  Cold  working  after 
oversize  drilling  would  be  an  example  of  increasing  a repair  cost 
and  obtaining  a better  repaired  crack  size  distribution. 

The  repair  costs  are  defined  in  the  input  as 
COST(J),  J=l,4.  Four  cost  values  are  required  even  if  some  are 
repeated . 

3.3  PROGRAM  PARAMETERS 

Two  other  parameters  must  be  specified  as  part  of  the  data 
input.  These  are  the  crack  size  increment,  DELTA,  and  the  number 
of  iterations  of  the  inspection  and  repair  cycle  to  be  performed, 

I CYCLE.  These  parameters  influence  only  the  operation  of  the 
program. 

DELTA  controls  the  number  of  points  at  which  the 
f“  (a)  and  f+(a)  functions  are  calculated.  Since  an  integration 
is  required  at  each  a value,  DELTA  should  not  be  chosen  arbitrarily 
small.  However,  computational  accuracy  can  be  lost  if  DELTA  is 
chosen  too  large.  A value  of  DELTA  approximately  equal  to  1/100 
of  the  crack  range  was  judged  to  provide  sufficient  accuracy  and 
acceptable  run  times. 

I CYCLE  merely  limits  the  number  of  iterations  of 
the  inspection/repair  calculations  to  the  desired  number.  The 
product  of  I CYCLE  and  UT  (inspection  interval)  determines  the 
fleet  life  being  covered  in  the  simulation. 


3.4  PROGRAM  LISTING 

The  program  users  manual  and  listing  are  presented  in 
Appendices  A and  B.  The  users  manual  is  on  a separate  card  deck 
and  does  not  necessarily  have  to  be  considered  as  a part  of  the 
program.  Appendix  A presents  instructions  for  calling  the  users 
manual  if  it  is  entered  in  the  computer,  a listing  of  the  manual, 
and  sample  output  from  one  cycle  of  a representative  computer  run. 

Appendix  B contains  the  complete  program  listing  for  PCRC 
(Predicted  Crack  Repair  Costs) . This  program  was  written  for  the 
CDC  6600  computer  located  in  the  Aeronautical  Systems  Division 
Computer  Center,  Area  B,  Wright-Patterson  Air  Force  Base.  It  is 
possible  that  differences  in  computer  systems  could  prevent  the 
direct  running  of  this  program  on  other  computers.  The  program 
also  utilizes  four  programs  from  the  International  Mathematical 
and  Statistical  Library  (IMSL) . These  routines  are  used  in  the 
calculation  of  integrals  ( DCADRE  and  DCSQDU)  and  for  inter- 
polation by  spline  function  fitting  (ICSICU)  and  evaluation 
(ICSEVU) . 

The  first  277  lines  of  code  control  the  input,  execution, 
and  output  of  the  program.  These  are  followed  by  the  definition 
of  8 subroutines,  of  which,  only  2 have  not  been  discussed  in  detail 
in  the  preceding  paragraphs.  Function  U2D  is  a subroutine  which 
evaluates  the  argument  in  the  calculation  of  FMINUS,  equation  (12). 
Function  FPLIN  controls  the  spline  fitting  used  in  interpolating 
to  obtain  intermediate  functional  values  in  the  numerical  integration 
routines . 

Results  from  the  programs  have  been  verified  under  a variety 
of  input  conditions.  However,  due  to  the  nature  of  the  numerical 
computations,  it  is  possible  that  inconsistent,  ill-defined,  or 
illogical  input  conditions  could  be  specified  which  would  cause 
failure  of  the  program.  Since  such  failures  would  not  necessarily 
preclude  completion  of  the  run,  the  output  should  always  be 
evaluated  to  insure  that  the  results  are  reasonable. 
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SECTION  IV 


EXAMPLE  RESULTS/SENSITIVITY  STUDY 

To  illustrate  the  use  of  the  Predicted  Crack  Repair  Cost 
computer  program,  a series  of  example  problems  were  run  which 
also  demonstrate  the  effect  of  variations  in  input  parameters  on 
the  predicted.  This  section  presents  an  example  for  defining 
an  input  data  set  and  the  results  of  the  limited  sensitivity 
study  which  was  conducted. 

The  approach  taken  in  the  sensitivity  study  was  to  define 
a set  of  "standard"  input  conditions  which  are  reasonably 
realistic  with  respect  to  Air  Force  experience  and  to  vary  only 
a small  portion  of  the  parameters  during  any  one  particular  set 
of  runs.  The  large  number  of  possible  combinations  of  the  input 
parameters  prohibited  a more  general  approach  in  which  all 
possible  combinations  are  studied.  Thus,  the  conclusions  of 
this  sensitivity  analysis  must  be  interpreted  in  this  restricted 
framework.  In  particular,  the  costs  of  inspection  and  repair 
can  greatly  influence  the  sensitivity  to  some  parameters  and  the 
allocation  of  a portion  of  the  total  cases  of  inspection  and 
repair  for  an  airplane  to  a single  structural  detail  is  not 
uniquely  defined.  If  a user  prefers  any  major  deviation  from 
the  input  "standard"  conditions,  a new  set  of  computer  runs  should 
be  generated  to  test  the  variables  of  interest. 

The  crack  growth  parameter  values  for  the  representative 
"standard"  set  were  derived  from  the  assumptions  that  the  critical 
crack  length,  aQ,  for  the  structural  details  was  0.24  in.,  that 
the  probability  of  an  equivalent  initial  crack  being  greater  than 
0.050  in. and  passing  the  manufacturers  quality  control  system 
was  0.001,  and  that  the  design  life  of  the  aircraft  was  6000 
flights.  Based  on  these  assumptions,  the  rationale  for  the 
assignment  of  specific  parameter  values  (in  the  order  requested 
by  the  program)  is  as  follows: 
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Inspection  interval:  UT=2000;  convenient  fraction  of 
design  life 

- Growth  rate:  UB=0 . 0002614 : from  equation  (34)  with  above 

assumptions 

- Usage  Bias:  UM=1.0;  design  usage 

Fit.  by  fit.  variability:  US=3.0;  representative  value 
observed  in  Miners  damage 

FQ  shape  parameter:  F0ALPn=2;  representative  values 
observed  in  F-4  data 

F0  scale  parameter:  FOBETA-O . 019 ; from  assumptions 
and  equation  (30) 

s^aPe  parameter:  FRALPH=2;  repaired  flow  distribution 
is  equivalent  to  initial 

FR  scale  parameter:  FRBETA=01 019 ; repaired  flaw  distri- 
bution is  equivalent  to  initial 
H shape  parameter:  HALPH=1;  representative  from 
inspection  reliability  data 

- H scale  parameter:  HBETA=0 . 1 ; assumed  value  for 

ultrasonic  roto-scanner  with  fastener  in 

- Crack  Length  Increment:  DELTA=0.004;  sufficient  resolution 

for  density  functions 

Number  of  Inspections:  ICYCLE=6;  perform  calculations 
to  12000  flight  lifetime 

Number  of  Items:  N=1000;  assume  a fleet  of  1000  aircraft 

- Inspection  flaw  sizes:AA(I) , 1=1,4  =0.000,  .016,  .032, 

.064;  repair  by  oversize  drilling  unless  crack 
is  greater  than  0.064  in.  which  then  requires 
replacement  of  entire  panel 

- Repair  Costs:  COST(j),  J=l,4=3,  3,  3,  100;  cost  of 

oversize  drilling  and  new  bolt  for  cracks  less 
than  0.064  in;  proportionate  share  of  panel 
replacement  (among  many  fasteners)  for  cracks 
over  0.064  in. 

Inspection  cost:  AIC=5000;  assume  proportionate  share 
of  ultrasonic  roto-scanner  inspection  with 
fastener  in  is  $5.00  per  detail  for  each  of  the 
1000  details. 
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Variations  from  the  above  set  of  standard  conditions 
were  investigated  for  inspection  intervals,  inspection  capabilities, 
repair  quality,  baseline  usage,  and  repair  costs.  In  addition, 
variations  in  equivalent  initial  and  repaired  flaw  quality  were 
analyzed  by  changing  shape  parameters  for  the  fixed  99.9 
percentile  and  also  by  reprogramming  two  other  families  of  distri- 
butions. The  results  of  these  studies  are  presented  in  the 
following  paragraphs.  To  display  the  changes  in  the  crack  size 
density  distributions,  Figure  5 displays  the  densities  for  the 
standard  conditions  immediately  after  the  repair  cycles. 

4.1  COMPARING  INSPECTION  INTERVALS 

To  compare  the  effect  of  time  between  inspection  intervals 
on  the  expected  costs  of  inspection  and  repair,  all  parameters 
were  fixed  at  standard  values  except  T.  In  separate  runs,  the 
cost  output  was  generated  for  inspection  intervals  of  1000,  2000, 

3000,  4000,  and  6000  flights.  The  results  are  summarized  in 

Table  II  which  presents  the  total  percentage  of  cracks  which  were 

identified  and  repaired,  the  expected  value  and  standard 

deviation  of  the  cost  of  each  inspection  and  repair  cycle,  and 

the  expected  value  and  standard  deviation  of  the  cumulative  inspection 

and  repair  costs.  The  calculations  were  carried  through  a 

lifetime  of  12000  flights  to  include  the  effects  of  utilization 

past  original  design  life.  Note  that  these  calculations  do  not 

include  any  expected  costs  which  could  arise  due  to  catastrophic 

failure  from  cracks  exceeding  the  0.24  in  critical  crack  size. 

The  results  of  Table  II  indicate  that  inspection  intervals 
of  2000,  3000,  4000  and  6000  flights  have  essentially  the  same 
expected  value  and  standard  deviation  of  total  costs  over  a 
lifetime  of  12000  flights  while  an  inspection  interval  of  1000 
flights  is  slightly  more  expensive.  Over  a lifetime  of  6000 
flights,  however,  there  appears  to  be  a cost  advantage  to  a 
3000  flight  inspection  interval.  Figure  6,  which  presents  the 
expected  total  costs  at  each  inspection  interval,  shows  these 
results  and  also  indicates  that  the  expected  cost  differences 

35 


- 


• .-.-r 


TABLE  II 


Variation  in  Time  Between  Inspections 
Standard  Values  for  Other  Parameters 


Inspection 

Interval 

2000 

3000 

Aircraft 

4000 

Life  (Flights) 

6000  8000 

9000 

10000 

12000 

1000 

P(i) 

0.219 

0.241 

0.254 

0.262 

0.260 

0.259 

0.258 

0.257 

Ci 

6200 

7600 

9200 

10900 

10900 

10800 

10700 

10700 

Si 

200 

400 

600 

700 

700 

700 

700 

700 

CTOT 

11800 

19400 

28600 

49900 

71700 

82500 

93200 

114700 

STOT 

200 

500 

800 

1200 

1600 

1700 

1900 

2100 

2000 

p<i> 

0.239 

6700 

0.  315 
14300 

0.360 

22000 

0.369 

23600 

0.362 

22600 

0.358 

22000 

si 

100 

900 

1200 

1200 

1200 

1200 

CTOT 

6700 

21000 

43000 

66600 

89200 

111200 

S 

TOT 

300 

900 

1500 

1900 

2300 

2600 

1000 

P(i) 

0.295 

0.425 

0.462 

0.444 

ei 

10800 

30600 

35800 

33000 

sjL 

CTOT 

700 

10800 

1400 

40400 

1500 

77200 

1400 

110200 

S 

TOT 

700 

1500 

2100 

2500 

4000 

fill 

0.361 

0.534 

0.526 

Ci 

19400 

46200 

43900 

si 

1100 

1600 

1500 

^TOT 

19400 

65600 

109500 

S 

TOT 

1100 

1900 

2 4 00 

6000 

0.515 

44200 

0.686 

66700 

Si 

1500 

1500 

CTOT 

44200 

110900 

S 

TOT 

1500 

2200 
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Inspection  Intervals 
Q - 1000  flights 
□ - 2000  flights 

- 3000  flights 

A ■ 4000  flights 
^ - 6000  flights 


o @ 


o o 

o □ 


between  intervals  of  2000,  3000,  4000,  and  6000  may  be  negligible. 


t 

! 
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The  higher  costs  for  the  1000  flight  inspection  interval 
can  be  attributed  to  the  relative  costs  of  the  inspection  to  the 
costs  of  the  repairs  being  encountered  at  this  interval.  Changing 
the  relative  magnitude  of  these  maintenance  costs  could  change 
the  conclusions  regarding  inspection  intervals.  For  example,  the 
high  cost  of  repairing  a crack  over  0.064  in.  in  length  is  due  to 
the  necessity  of  replacing  an  entire  panel.  If  the  cost  of  the 
panel  is  $20,000  and  the  entire  cost  must  be  allocated  to  the 
critical  location  being  simulated,  then  C (4) =20, 000  and  the 
expected  values  and  standard  deviations  for  various  inspection 
intervals  up  to  a 6000  flight  life  would  appear  as  in  Table  III. 
These  results  clearly  indicate  the  significant  advantage  of  the 
shorter  inspection  intervals  over  the  longer  as  well  as  the 
greatly  increased  expected  maintenance  costs.  Note,  however, 
that  the  data  of  Table  II  pertain  to  only  one  fastener  hole  while 
the  panel  replacement  which  causes  the  high  costs  of  Table  III 
repairs  several  hundred  fastener  holes. 

4.2  COMPARING  INSPECTION  CAPABILITIES 

To  compare  inspection  capabilities,  three  types  of 
inspections  were  hypothesized  with  different  relative  costs  and 
•relative  capabilities.  These  are  summarized  as  follows: 


Inspection  Type 


Relative  Cost/ 
Detail 


Relative  Capability 
(Characteristic  Value) 


Visual-fastener  in  $1.25 

Roto-scanner  Ultrasonic-  $5.00 

fastener  in 

Motorized  eddy  current-  $12.50 

fastener  out 


0.2  in. 
0.1  in . 

0.05  in. 


These  values  were  judged  to  be  representative  of  Air  Force  experience. 
(The  roto-scanner  ultrasonic  parameter  values  are  "standard"  as 
previously  defined) . 


TABLE  III.  Variation  in  Time  Between  Inspections  - High 
Cost  for  Repair  of  Cracks  over  0.064  in. 


Inspection 

Aircraft  Life  (Flights) 

Interval 

/ 000 

2000 

3000 

4000 

5000 

6000 

1000 

Ci 

18,000 

119,500 

385,600 

714,500 

952,200 

1,056,400 

Si 

15,700 

47,600 

86,300 

116,900 

134,300 

141,100 

c~" 

TOT 

18,000 

137,500 

523,100 

1,237,600 

2,189,800 

3 ,246,200 

STOT 

15,700 

50,100 

99,800 

153,800 

204,100 

248,200 

2000 

r* 

“i 

199,300 

1,720,000 

3,288,000 

Si 

61,900 

177,100 

234,300 

CTOT 

199,300 

1,920,000 

5,203,000 

STOT 

61,900 

187,600 

300,100 

3000 

~i 

1,009,300 

5,011,100 

Si 

138,100 

274,000| 

CTOT 

1,009,300 

6,020,400; 

STOT 

138,100 

306,300 

6000 

=7 

7,777,600 

Si 

308, 3001 

CTOT 

7,777,600! 

STOT 

308,300 

The  results  using  the  costs  of  repair  from  the  standard 
set  of  parameters  are  presented  in  Table  IV.  For  these  relative 
costs  of  inspections  and  repairs,  the  expected  maintenance  costs 
are  less  for  the  cheaper,  less  effective  inspection  procedures. 

This  cost  savings  is  due  to  the  reduced  cost  of  the  inspections 
and  may  also  reflect  a false  savings  of  not  detecting  (and  not 
repairing)  the  cracks  greater  than  0.064  in.  Note  that  these 
cost  figures  do  not  reflect  expected  losses  that  may  occur  when 
the  cracks  reach  the  critical  size  of  0.24  in.  Due  to  the 
assumptions  which  determine  the  standard  parameters,  less  than 
0.1  percent  of  the  aircraft  will  have  a crack  size  greater  than 
0.24  in.  at  6000  flights. 

Table  V presents  the  expected  cost  data  when  the  repair 
cost  for  a crack  longer  than  0.064  in.  is  $20,000.  At  the  6000 
flight  aircraft  life,  the  best  and  highest  cost  inspection 
procedure  is  the  most  cost  effective.  At  the  earlier  inspection, 
the  high  cost  of  the  expensive  inspection  procedure  is  attributable 
to  the  discovery  of  more  of  the  large  cracks.  Over  the  6000 
flight  life  of  the  aircraft,  the  early  repair  of  these  cracks 
becomes  evident.  Again,  the  addition  of  expected  loss  due  to 
catastrophic  failures  would  reinforce  the  cost  efficiency  of 
the  expensive  inspection  procedure. 

4.3  COMPARING  REPAIR  QUALITY 

To  compare  repair  quality,  it  was  assumed  that  after 
oversizing  the  fastener  hole,  it  was  either  cold  worked  or 
a fatigue  fastener  was  installed.  It  was  assumed  that  this 
operation  would  cost  twice  as  much  ($6  per  repair)  for  cracks 
less  than  0.064  in.  and  that  the  equivalent  repair  quality 
distribution  would  have  a characteristic  value  of  0.0019  (ten 
times  better).  Repair  of  cracks  greater  than  0.064  in.  are 
accomplished  by  panel  replacement  and  would  cost  the  same. 

Table  VI  presents  the  summary  of  the  maintenance  cost  data 
for  the  standard  cost  of  repair  of  cracks  greater  than  0.064  in. 
(i.e.,$100).  In  a 6000  flight  lifetime,  the  more  expensive  repair 
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TABLE  IV 


Variation  in  Inspection  Capabilities 
Standard  Values  for  Other  Parameters 


Inspection 

Capabilities 

2000 

Aircraft  Life  (Flights) 

4000  6000  8000 

10000 

12000 

HBETA-0. 2 

P(i) 

0.130 

0.190 

0.249 

0.285 

0.289 

0.227 

I(c)»1250 

ST 

2200 

8100 

17000 

21900 

22000 

20300 

S. 

1 

240 

780 

1100 

1300 

1300 

1200 

CT0T 

2200 

10300 

27300 

49200 

71200 

91500 

ST0T 

240 

820 

1400 

1900 

2300 

2600 

HBETA-0.1 

P(i) 

0.239 

0.315 

0.360 

0.369 

0.362 

0.357 

I(C)»5000 

6700 

14300 

22000 

23600 

22600 

22000 

si 

300 

900 

1200 

1200 

1200 

1200 

CTOT 

6700 

21000 

43000 

66600 

89200 

111200 

STOT 

300 

900 

1500 

1900 

2300 

2600 

HBETA-0. 05 

P (i) 

0.408 

0.471 

0.486 

0.484 

0.481 

0.480 

1(0-12500 

15100 

21700 

25100 

25100 

24900 

24800 

si 

380 

860 

1000 

1000 

1000 

1000 

CT0T 

15100 

36800 

61900 

87000 

111900 

136700 

s 

TOT 

380 

940 

1400 

1700 

2000 

2200 
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TABLE  V.  Variation  in  Inspection  Capabilities- 
High  Cost  for  Repair  of  Cracks  Over 
0.064  in. 


Inspection 

Capabilities 

Aircraft 

2000 

Life  (Flights) 

4000 

6000 

HBETA-0. 2 

Ci 

115,600 

1,302,600 

3,097,300 

1(0-1250 

s. 

1 

47,600 

156,000 

228,800 

CTOT 

115,600 

1,431,200 

4,515,500 

ST0T 

47,600 

163,100 

280,900 

HBETA-0. 1 

199,300 

1,720,200 

3,288,000 

1(0-5000 

Si 

61,900 

177,100 

234,300 

CT0T 

199,300 

1,920,000 

5,208,000 

STOT 

61,900 

187,600 

300,100 

HBETA-0. 05 

302,000 

1,618,800 

2,317,700 

I (0-12500 

s. 

1 

75,400 

171,800 

201,900 

CT0T 

302,000 

1,920,800 

4,238,500 

STOT 

75,400 

187,600 

275,600 

TABLE  VI.  Variations  in  Repair  Quality  Standard 
Values  for  Other  Parameters 


Repair 

Quality 

2000 

Aircraft  Life 

4000 

(Flights) 

6000 

8000 

10000 

12000 

RBETA-0.019 

P(i> 

0.239 

0.315 

0.360 

0.369 

0.362 

0.358 

C(J)-3 

^7 

6700 

14300 

22000 

23600 

22600 

22000 

S. 

1 

300 

900 

1200 

1200 

1200 

1200 

c 

TOT 

6700 

21000 

43000 

66600 

89200 

111200 

STOT 

300 

900 

1500 

1900 

2300 

2600 

RBETA-0.0019 

P (i) 

0.239 

0.260 

0.228 

0.156 

0.088 

0.047 

C(J)-6 

7300 

14400 

19600 

17000 

11600 

7900 

Si 

300 

900 

1100 

1000 

300 

500 

CTOT 

7300 

21700 

41300 

58300 

69900 

77800 

STOT 

300 

900 

1400 

1800 

1900 

2000 

is  only  slightly  less  expensive  but  beyond  6000  flights  the 
effect  of  the  improved  quality  of  repair  becomes  significant. 

Table  VII  summarizes  the  repair  cost  data  when  the  repair  cost 
(panel  replacement)  of  the  cracks  longer  than  0.064  in.  is 
$20,000,  i.e.,the  total  cost  is  attributed  to  the  critical 
location.  The  cost  effectiveness  of  the  higher  cost,  higher 
quality  repair  under  these  conditions,  is  significant  at  the 
6000  flight  lifetime. 

4.4  COMPARING  DIFFERENCES  IN  AIRCRAFT  USAGE 

Differences  in  aircraft  usage  are  modeled  by  the  usage 
bias  and  flight  by  flight  variability  parameters.  To  determine 
the  sensitivity  of  predicted  maintenance  costs  to  these  parameters, 
the  program  was  exercised  for  usage  bias  values  of  0.3,  0.9,  1.0, 
1.1,  and  1.2  and  for  flight  by  flight  variability  (coefficients 
of  variation)  of  1.5,  3.0  and  4.5.  The  maintenance  costs  were 
calculated  for  both  the  high  and  the  low  panel  replacement  costs 
(repair  costs  for  cracks  larger  than  0.064  in). 

The  effect  due  to  usage  bias  is  signifcant  as  expected. 

Table  VIII  shows  the  expected  total  inspection  and  repair  costs  at 
each  inspection  for  both  sets  of  panel  replacement  costs.  The 
effect  of  deviations  from  the  baseline  are  proportionately  larger 
for  the  higher  repair  costs.  If  a baseline  usage  is  assumed  when 
the  aircraft  are  being  used  in  a 10  percent  more  severe  environment, 
the  expected  maintenance  costs  over  a 6000  flight  life  will  be 
underestimated  by  26  percent  for  the  high  repair  costs  and  by 
12  percent  for  the  low  repair  cost. 

The  effect  of  flight  by  flight  variability  on  the  expected 
total  costs  was  negligible.  The  range  of  values  considered 
(1.5  to  4.5)  encompasses  the  variability  that  would  be  expected 
in  operation  and  the  expected  costs  were  within  1 percent  of 
baseline  for  both  sets  of  cost  of  repair.  This  result  agrees  with 
a similar  conclusion  regarding  flight  by  flight  variability 
of  Miners  damage  (Reference  12).  However,  even  though  flight  by 
flight  variability  may  be  unimportant  in  the  calculation  of 
expected  costs,  it  is  a necessary  parameter  for  the  estimation 
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TABLE  VII.  Variation  in  Repair  Quality-High  Cost 
for  Repair  of  Cracks  Over  0.064  in. 


Repair 

Quality 

Aircraft 

2000 

Life  (Flights) 

4000 

6000 

RBETA-0.019 

Ci 

199,300 

1,720,700 

3,238,000 

C(J)-3 

Si 

61,900 

177,100 

234,000 

CTOT 

199,300 

1,920,000 

5,208,000 

STOT 

61,900 

187,600 

300,100 

RBETA*0. 0019 

200,000 

1,674,800 

2,828,700 

C ( J)  =6 

si 

61,900 

174,900 

220,200 

c 

TOT 

200,000 

1,874,800 

4,703,500 

STOT 

61,900 

185,500 

283,000 

TABLE  VIII.  Variation  in  Usage  Bias-Expected 
Total  Maintenance  Cost 


Panel 

Aircraft  Life 

(Flights) 

Replacement 

Usage 

Cost 

Bias 

2000 

4000 

6000 

C (4) “$100 

0.8 

6000 

15700 

30900 

0.9 

6300 

18000 

36600 

1.0 

6700 

20900 

43000 

1.1 

7100 

24300 

49600 

1.2 

7800 

28100 

56500 

C(4)-$20,000 

mm 

82,600 

373,700 

2,791,500 

l 

130,600 

1,340,000 

3,938,300 

1 

199,300 

1,920,000 

5,208,000 

293,500 

2,601,700 

6,554,300 

i 

K£ 1 

417,800 

3,368,600 

7,939,300 

of  the  precision  of  the  usage  bias  estimate.  Thus,  the  importance 
of  this  parameter  is  in  the  determination  and/or  the  evaluation 
of  sample  size  considerations  with  respect  to  the  estimate  of 
the  significant  parameter  of  usage  bias. 

4.5  COMPARING  FLAW  SIZE  DISTRIBUTIONS 

The  equivalent  flaw  size  distributions  for  either  new  or 
repaired  structure  cannot  be  ascertained  by  direct  observations 
of  cracks.  Rather,  observations  representative  of  this  distri- 
bution are  obtained  from  a back  calculation  of  an  observed 
crack  to  time  zero.  The  crack  may  or  may  not  have  been  present 
at  time  zero  but  the  structure  responded  as  though  it  were  and, 
hence,  the  equivalence  of  the  crack  at  time  zero. 

Due  to  the  difficulty  in  obtaining  valid  equivalent 
time  zero  crack  lengths,  relatively  few  such  data  points  have 
been  generated.  There  is  very  little  data  available  for  use 
in  selecting  an  appropriate  family  of  distributions  for 
describing  initial  flaw  sizes.  As  one  step  in  assessing  the 
relative  importance  of  the  possible  shapes  of  this  distribution, 
expected  inspection  and  repair  costs  were  calculated  for  three 
different  families  of  distributions  with  three  different  medians 
(and,  hence,  shapes)  for  each  family. 

The  three  families  of  distributions  that  were  compared  for 
initial  and  repaired  equivalent  flaw  size  distributions  were 
the  Weibull,  the  gamma,  and  the  lognormal.  The  equation  of  the 
Weibull  density  function  is  given  in  equation  (29) . The 
equation  of  the  density  function  of  the  gamma  distribution  is 
given  by 


f (y) 


a a a-1 

I 1 

T (a) 


exp 


[6  y] 


while  that  of  the  lognormal  density  is  given  by 


(35) 


f (y) 


1 

y ol/2t\ 


exp 


(in  y-f3)2 

2a2 


(36) 
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For  each  of  the  three  families  the  parameters  were  determined  such 
that  99.9  percent  of  the  flaws  were  less  than  0.050  in.  and  they 
had  equal  medians  (50th  percentiles) . The  parameter  values  used 
are  given  in  Table  IX.  The  resulting  density  functions  are 
pictured  in  Figure  7.  The  distributions  with  a median  of  0.0156 
have  a shape  parameter  approximately  equal  to  that  of  the  A-7D 
initial  quality  data.  Reference  8. 

The  expected  total  costs  of  inspection  and  repair  are 
presented  in  Table  x for  each  of  the  three  families,  with  the 
three  different  median  initial  (and  repaired)  cracks,  and  for  the 
high  and  low  cost  for  repair  of  cracks  longer  than  0.064  in. 

The  differences  in  the  expected  costs  can  be  interpreted  with 
the  aid  of  the  density  functions  of  Figure  7.  The  Weibull 
density  function  for  the  three  medians  (or  the  three  shapes) 
considered  is  more  diffuse  and,  hence,  has  greater  area  in  the 
tails  than  the  other  two  distributional  types.  Similarly,  the 
gamma  distribution  is  more  diffuse  than  the  lognormal.  In  the 
early  life  of  the  aircraft,  these  larger  cracks  will  grow  more 
rapidly,  have  a greater  chance  of  being  detected,  and  hence, 

will  incur  higher  repair  costs.  However,  when  cracks  around  the 
central  tendencies  of  the  distributions  have  grown  to  more 
detectable  sizes,  the  other  distributions  will  have  a larger 
percentage  requiring  repair.  Hence,  the  expected  costs  cross 
over  at  an  inspection  time  which  depends  on  the  median.  The 
larger  the  median,  the  sooner  the  crossover  occurs. 

Note  that  for  the  lower  cost  of  repair  of  the  long 
cracks,  the  differences  in  expected  total  costs  for  the  three 
families  is  not  practically  significant.  Further,  the  differences 
for  the  higher  costs  of  repair  are  not  excessively  large  at 
6000  flights  for  the  two  higher  medians  or  for  any  of  the  medians 
at 8000  and  10000  flight  lifetimes.  However,  the  differences  in 
expected  total  costs  for  the  three  medians  is  significant  for 
all  three  families  at  each  inspection  throughout  aircraft  life. 
Therefore,  from  the  viewpoint  of  maintenance  costs  it  would  be 
more  important  to  determine  the  shape  (or  median)  of  the  density 
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TABLE  IX.  Parameter  Values  for  Comparing  Families 
of  Initial  Flaw  Size  Distributions 


Median 

Weibull  gamma  lognormal 

0.0091  a 

e 

1.347  2 0.552 

0.11S1  185  -4.703 

0.0156  a 

8 

1.971  5 0.378 

0.01875  300  -4.162 

0.0210  a 

8 

2.653  10  0.280 

0.02413  460  -3.862 

Log  Normal 


Figure".  Equivalent  Initial  Crack  Size  Density  Functions. 


TABLE  X. 


Comparing  Families  of  Equivalent  Flaw  Size 
Distr ibutions-Expected  Total  Maintenance  Costs 


Pane  1 

Flaw 

Aircraft  Life  (Flights) 

Replacement 

Size 

Cost 

Median 

Distribut ion 

2000 

4000 

6000 

8000 

10000 

C ( 4 ) *$ 100 

0.0091 

W 

5900 

14700 

28300 

45000 

62300 

G 

5800 

14200 

27300 

44300 

62300 

LN 

5700 

13100 

25000 

42800 

62800 

0.0156 

w 

6600 

20600 

42300 

65700 

88100 

G 

6100 

18900 

41200 

66200 

89400 

LN 

6200 

17800 

40800 

67000 

90400 

0.0210 

w 

7700 

28700 

54700 

85000 

112500 

G 

7000 

27100 

58100 

87500 

1 14900 

LN 

6800 

26500 

58500 

88100 

115600 

C ( 4 ) =$20,000 

0.0091 

w 

96000 

1,865,000 

2,  370,000 

4,610,000 

6,973.000 

G 

81500 

1,502,000 

2,166,000 

4 , 457 , 000 

6 , 955 , 000 

LN 

57000 

1,260,000 

1,668,000 

4,121 ,000 

7,015,000 

0.0156 

w 

195,000 

3,450,000 

5,087,000 

R , 657 , 000 

12.01 1,000 

G 

128 , 000 

3,101,000 

4,849,000 

8,743,000 

12,274,000 

LN 

104,000 

2,989,000 

4,746,000 

8,883,000 

12,438, 000 

0.0210 

W 

381,000 

3,450,000 

8,076, 000 

12,606,000 

1 6 , 117  7 , 000 

G 

222,000 

3,101,000 

8,184, 000 

12, 94  ‘>,000 

17,111, 000 

LN 

190,000 

2,989,000 

8,275,000 

1 1,070,000 

l 1 , 4 4H , 000 

function  of  equivalent  flow  sizes  than  to  determine  the  family 
of  distributions.  Note,  however,  that  safety  of  flight  considera- 
tions which  are  based  only  on  the  largest  equivalent  initial 
flaw  sizes  may  well  indicate  the  reverse  order  of  priority. 


SECTION  V 


SUMMARY  AND  RECOMMENDATIONS 

There  were  three  primary  objectives  of  this  study:  (1)  to 
write  a computer  program  which  could  be  used  to  predict  the  repair 
costs  of  cracks  which  develop  during  operational  aircraft  usage, 

(2)  to  provide  a document  which  describes  the  use  of  the  computer 
program  and  guides  a potential  user  in  the  specification  of 
input  requirements,  and  (3)  to  use  data  which  is  reasonably  repre- 
sentative of  Air  Force  experience  as  input  for  use  in  determining 
the  sensitivity  of  predicted  crack  repair  costs  to  variations  in 
the  input  parameters. 

The  statistical  model  which  has  been  programmed  to  achieve  these 
objectives  computes  at  each  inspection  the  crack  size  densities 
before  the  inspection  and  after  the  repair  of  cracks  which  have 
been  found;  the  percentage  of  cracks  in  four  ranges  of  crack  size; 
the  expected  value  and  standard  deviation  of  the  cost  of  each 
inspection  and  repair  cycle;  and  the  expected  value  and  standard 
deviation  of  the  total  costs  of  all  inspection  and  repair  cycles. 

The  crack  size  distribution  output  is  optional.  The  model  is  pro- 
grammed with  a Weibull  distribution  for  the  equivalent  initial  and 
repaired  crack  sizes  and  for  the  probability  of  crack  detection  as 
a function  of  crack  size.  It  is  assumed  that  cracks  grow  in  the 
operational  environment  in  accordance  with  an  exponential  model  as 
derived  from  a baseline  spectrum,  and  that  deviations  from  baseline 
crack  growth  and  flight  by  flight  usage  variability  can  be 
accounted  for  by  a normal  distribution.  The  model  allows  for  different 
repair  costs  for  up  to  four  crack  size  ranges  with  the  largest 
range  usually  accounting  for  the  cost  of  panel  replacement  or, 
possibly,  the  cost  of  catastrophic  failure.  Changes  in  the  family 
of  distributions  for  the  crack  sizes  and  the  probability  of  crack 
detection  can  be  easily  accomplished  if  the  replacement  family  is 
defined  in  terms  of  1 or  2 parameters. 
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The  program  was  demonstrated  by  computed  expected  maintenance 
costs  for  different  inspection  intervals,  inspection  capabilities, 
repair  quality,  aircraft  usage,  and  equivalent  flaw  size  distribu- 
tions. The  computations  were  made  for  a set  of  input  data  that  is 
reasonably  representative  of  Air  Force  maintenance  capabilities 
and  costs.  However,  the  risk  of  catastrophic  failure  was  not 
included.  The  conclusions  of  the  sensitivity  studies  are  dependent 
on  the  values  of  the  other  parameters  and  only  relatively  few 
parameter  combinations  could  be  tested.  The  sensitivity  computa- 
tions were  performed  for  a high  and  low  cost  of  repairing  long  cracks 
since  most  results  were  definitely  sensitive  to  this  parameter. 

The  following  general  conclusions  resulted: 

1)  The  effect  of  inspection  interval  depended  on  the  cost 
of  repairing  long  cracks.  For  the  high  repair  cost,  frequent 
inspections  saved  money  but  for  the  low  repair  cost  of  long  cracks, 
a tradeoff  between  cost  of  inspection  and  cost  of  repair  resulted 
in  cost  savings  at  the  less  frequent  inspection  intervals. 

2)  The  effect  of  inspection  capabilities  is  highly  dependent 
on  the  relative  magnitudes  of  the  cost  of  inspection  and  cost  of 
repair  of  long  cracks.  For  the  cost  values  used  in  this  study, 
cheaper,  less  effective  inspections  resulted  in  significant 

(35  percent)  cost  savings  for  the  low  repair  cost.  For  the  high 
repair  cost,  however,  the  most  expensive  (most  effective)  inspection 
led  to  the  least  total  maintenance  cost  while  the  intermediate 
inspection  led  to  the  highest  total  maintenance  cost. 

3)  The  two  repair  quality  capabilities  tested  indicated 
that  the  more  expensive,  more  effective  repair  is  cost  effective 
at  the  design  life. 

4)  Deviations  of  10  percent  or  more  from  baseline  usage  are 
significant  with  differences  in  expected  maintenance  costs  of  about 
15  percent  being  observed  for  a 10  percent  change  in  usage  severity. 
Flight  by  flight  variability  in  usage  severity,  however,  has  no 
practically  significant  effect  on  expected  maintenance  costs 

(less  than  1 percent) . 


5)  Three  families  of  distributions  for  equivalent  initial 
flaw  sizes  were  tested.  The  differences  in  expected  total  costs 
at  the  design  life  from  the  three  families  were  negligible  in 
comparison  to  the  location  of  the  median  flaw  size  within  any  one 
of  the  families.  Therefore,  it  is  more  important  to  properly 
locate  the  distribution  than  to  be  concerned  about  the  statistical 
mode 1 . 

With  respect  to  recommendations,  there  are  two  directions  for 
further  work  that  are  indicated  as  a result  of  this  study.  First, 
the  program  can  be  modified  to  permit  input  being  defined  in  more 
qeneral  terms.  For  example,  it  may  be  desirable  to  specify  different 
parameter  values  at  every  inspection  interval  for  every  input  item. 
Further,  it  may  be  desirable  to  pursue  a random  growth  concept  by 
allowing  the  crack  growth  parameters  of  the  model  to  be  a random 
variable.  This  concept  would  require  changing  the  probability 
calculations  of  the  model  and  would  also  necessitate  some  investiga- 
tion as  to  the  type  of  distribution  that  would  be  used  to  model  this 
parameter . 

The  second  direction  suggested  by  the  study  is  to  pursue 
the  usage  parameter  as  a method  for  monitoring  for  change  in  usage. 
This  direction  would  require  the  processing  of  stress  histories 
(perhaps,  from  different  aircraft  types)  and  the  calculation  of 
c rack  length  as  a function  of  usage  time.  This  direction  is  a 
distinct  change  from  the  current  study  but  could  have  significance 
in  other  force  management  applications. 


APPENDIX  A 
USERS  MANUAL 


IN  ORDER  TO  LIST  THE  USER  MANUAL , AFTER  A COMMAND  REQUEST 
THE  USER  SHOULD  ENTER? 

ATTACH , FILE1 , PCRCDOCUMENTATION 
EDITOR 

THE  USER  IS  NOW  IN  THE  EDITOR  MODE*  AFTER  A DOUBLE  DOT  IS 
RETURNED » ENTER? 


EDIT,FILE1,S 
LIST » ALL 


THE  LISTING  WILL  LOOK  LIKE  THIS? 


USER  MANUAL 

(NOTE?  REFERBACK  FOLLOWED  BY  A NUMBER  (1  THROUGH  9)  IS  FOR 
USER  REFERENCE  WITHIN  THE  USER  MANUAL  AND  SHOULD  NOT  BE 
MISCONSTRUED  AS  PART  OF  THE  INPUT  OR  OUTPUT  OF  THE  PROGRAM.) 

AFTER  LOGIN  PROCEDURES,  THE  USER  MUST  INITIATE  EXECUTION  OF 
THE  PROGRAM*  TO  DO  THIS,  THE  USER  MUST  INPUT  EACH  OF  THE 
FOLLOWING  WHEN  A COMMAND  IS  REQUESTED? 

ATTACH, FILE2,PCRC 

ATTACH,LIB, IMSL, ID=LIBRARY,SN=ASD 

LIBRARY, LIB 

RETURN  OUTPUT 

CONNECT  TAPE9 

FTN, I=FILE  2,  L=0 

OUTPUT  IS  RETURNED  (DISCONNECTED)  TO  BYPASS  PRINTING  OF 
NUMEROUS  TRIVIAL  ERROR  WARNINGS.  AFTER  PROCESSING  OF  PCRC 
IS  COMPLETED,  THE  USER  MUST  ENTER?  CONNECT  OUTPUT 

AFTER  A COMMAND  REQUEST  IF  PROCESSING  IS  TO  PROCEED  ON  OTHER 
PROGRAMS. 


IF  THE  USER  ALSO  WANTS  TO  SEE  THE  CRACK  SIZE  DENSITY  FUNCTIONS 
AS  PART  OF  THE  OUTPUT,  THE  USER  WILL  ADDITIONALLY  INPUT  AFTER  A 
COMMAND  REQUEST? 

CONNECT  TAPE! 

THE  PROGRAM  WILL  BEGIN  EXECUTION  AT  THIS  POINT  WHEN  THE  USER 
ENTERS  AFTER  A COMMAND  REQUEST?  LGO 

THIS  PAGE  IS  BEST  QUALITY  PRACtlQAfiUI 
KftUM  00 1 Y n/KAlSHED  TO  DOC  „ 


in-  > ihiit  - 


Tr!E  USER  MUST  BE  VERY  CAREFUL  IN  ENTERING  RESPONSES. 

FOR  INSTANCE#  THE  USER  MUST  ENTER  YESr  Yr  NO#  N#  OR  VARIABLE 
NAMES  WITHOUT  ANY  BLANKS  TO  INSURE  ACCURATE  READING  OF 
THE  RESPONSE  BY  THE  PROGRAM.  THE  PROGRAM  BEGINS 
BY  OUTPUTTING  THE  FOLLOWING. 

(REFERBACK  1)J 

ENTER  UT»UB#UM#US. . i 

THE  USER  ENTERS  VALUES  FOR  EACH  FIELD  SEPARATED  BY  COMMAS# 

WHERE:  UT  IS  INSPECTION  INTERVAL# 

UB  IS  GROWTH  RATE# 

UM  IS  USAGE  BIAS#  AND 

US  IS  FLT.  BY  FLT . VARIABILITY. 

THE  PROGRAM  RESPONDS  WITH: 

ENTER  FOALPH#  FOBETA..: 

THE  USER  ENTERS  VALUES  FOR  EACH  FIELD  SEPARATED  BY  COMMAS# 

WHERE?  FOALPH  IS  SHAPE  PARAMETER  AND 

FOBETA  IS  SCALE  PARAMETER 
FOR  THE  INITIAL  FLAW  SIZE  DISTRIBUTION. 

THE  PROGRAM  RESPONDS  WITH? 

ENTER  FRALPH#FRBETA. . : 

t 

THE  USER  ENTERS  VALUES  FOR  EACH  FIELD  SEPARATED  BY  COMMAS# 

WHERE:  FRALPH  IS  SHAPE  PARAMETER  AND 

FRBETA  IS  SCALE  PARAMETER 
FOR  THE  REPAIRED  FLAW  SIZE  DISTRIBUTION. 

THE  PROGRAM  RESPONDS  WITH: 

ENTER  HALPH  r HBETA . . : 

THE  USER  ENTERS  VALUES  FOR  EACH  FIELD  SEPARATED  BY  COMMAS# 

WHERE.*  HALPH  IS  SHAPE  PARAMETER  AND 
HBETA  IS  SCALE  PARAMETER 
FOR  THE  PROBABILITY  OF  CRACK  DETECTION. 

THE  PROGRAM  RESPONDS  WITH: 

Em  TER  DELTA# I CYCLE » N. . : 

THE  USER  ENTERS  VALUES  FOR  EACH  FIELD  SEPARATED  BY  COMMAS# 
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where:  delta  is  increment  for  crack  size  density  functions 

ICYCLE  IS  NUMBER  OF  INSPECTIONS  AND 
N IS  NUMBER  OF  ITEMS, 

THE  PROGRAM  RESPONDS  WITH: 


(REFEREACK  2): 

ENTER  AA  < I ) * 1 = 1 > 4 . « * 

THE  USER  ENTERS  THE  VALUES  FOR  THE  FOUR  ELEMENTS  OF  THE  ARRAY 
WHERE*.  AA(I)  DEFINE  LEVEL  I CRACKS 

THE  PROGRAM  RESPONDS  WITH: 

< REFERBACK  3): 

ENTER  COST ( J ) » J=1 v A . . l 

THE  USER  ENTERS  THE  VALUES  FOR  THE  FOUR  ELEMENTS  OF  THE  ARRAY 
WHERE:  COST(J)  IS  THE  REPAIR  COST  FOR  LEVEL  J CRACKS, 

THE  PROGRAM  RESPONDS  WITH4, 


(REFERBACK  4 > 4. 

DO  YOU  WANT  MULTIPLE  VALUES  FOR  AIC..‘. 

IF  THE  INSPECTION  COSTS  WILL  VARY  FOR  EACH  ICYCLE 
THE  USER  ENTERS  YES  OR  Y.  THE  PROGRAM  RESPONDS  WITH: 

ENTER  VALUES  FOR  AIC— ONE  VALUE  PER  ICYCLE.,: 

THE  USER  ENTERS  AS  MANY  VALUES  FOR  AIC  AS  THE  VALUE  OF  ICYCLE. 
IF  THE  INSPECTION  COSTS  WILL  REMAIN  THE  SAME  FOR  EACH  ICYCLE 
THE  USER  ENTERS  NO  OR  N.  THE  PROGRAM  RESPONDS  WITH,4 

ENTER  VALUE  FOR  AIC.,: 

THE  USER  ENTERS  ONE  VALUE » THE  PROGRAM  USES  THIS  VALUE  FOR 
EACH  ICYCLE. 


( REFERBACK  5): 


TR  . 
TH  . 
Ti  ; 


PROGRAM  WILL  RESPOND  WITH  A LISTING 
:•  VARIABLE  NAMES*  AND  WHAT  THE  VAR! 
IVOGFAM  WILL  THEN  ASIC 


OF  ALL  INPUT  VALUES 
ABLE  NAME  REPRESENTS 


» 
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IS  CRACK  SIZE  DATA  CORRECT . « J 


THE  USER  ENTERS  YES  OR  Y IF  CALCULATIONS  ARE  TO  PROCEED. 

IF  THE  USER  WANTS  TO  CHANGE  SOME  INPUT  DATA,  THE  USER  ENTERS 
NO  OR  N.  THE  PROGRAM  WILL  BYPASS  EXECUTION  OF  CALCULATIONS 
FuR  THIS  DATA  AND  ALLOWS  THE  USER  TO  CHANGE  DATA  AS  IF  THE 
USER  HAS  REQUESTED  ADDITIONAL  PROCESSING  (REFERBACK  7). 

ASSUMING  THE  USER  WANTS  CALCULATIONS  STEPS  TO  BE  EXECUTED, 

THE  PROGRAM  WILL  PRINT  OUT  A SET  OF  CALCULATIONS  FOR  EACH 
ICYCLE.  IF,  IN  ADDITION,  THE  USER  HAS  REQUESTED  THE  DENSITY 
FUNCTION,  ITS  VALUES  WILL  BE  DISPLAYED  BETWEEN  SETS  OF 
CALCULATIONS.  AT  THE  END  OF  THE  SET  OF  CALCULATIONS,  THE  PROGRAM 
WILL  ask: 


(REFERBACK  6)1 

DO  YOU  WANT  TO  CONTINUE  PROCESSING..: 

IF  NO  OR  N IS  ENTERED  THE  PROCESSING  OF  PCRC  STOPS. 

IF  YES  OR  Y IS  ENTERED  OR  IF  THE  CALCULATIONS  HAVE  BEEN 
BYPASSED  THE  PROGRAM  ASKS: 

(REFERBACK  7)1 


DO  YOU  WANT  TO  RE-INITI ALIZE  ALL  VALUES,.: 

IF  YES  OR  Y IS  ENTERED,  THE  PROGRAM  WILL  ASK  FDR  INPUT  DATA 
LINE  BY  LINE  AS  IN  THE  INITIATION  OF  THE  PROGRAM  (REFERBACK  1) 
AND  BYPASS  FURTHER  QUESTIONS  FOR  RE-INITIATION. 

IF  NO  OR  N IS  ENTERED,  THE  PROGRAM  WILL  ASK: 

DO  YOU  WANT  TO  RE-INITI ALIZE  THE  VALUES  OF  AA<I).\  : 

IF  YES  OR  Y IS  ENTERED,  THE  PROGRAM  WILL  ASK  FOR  VALUES  AS  IN 
THE  INITIAL  READ  FOR  AA(I>  (REFERBACK  2). 

REGARDLESS  OF  THE  PRECEEDING  ANSWER,  THE  PROGRAM  WILL  ASKJ 

DO  YOU  WANT  TO  RE- I MI T I ALIZE  THE  VALUES  OF  COST(J)..: 


IF  YES  OR  Y IS  ENTERED,  THE  PROGRAM  WILL  ASK  FOR  VALUES  AS  IN 
THE  INTIAL  READ  FOR  COST(J)  (REFERBACK  3). 

REGARDLESS  OF  THE  PRECEEDING  ANSWER,  THE  PROGRAM  WILL  ASKJ 


ENTER  THE.  VARIABLE  NAME  WHOSE  VALUE  IS  TO  BE  RE-INITIALIZED.. 


I ;• . 


..  rt  if .. 


SHOULD  BE  AWARE  THAT  THE 
NAME  ENTERED.  THE  USER 
i-.'Ar.ZS  EXACTLY  AS  LISTED: 


PROGRAM  SEARCHES  FOR  THE 
SHOULD  ENiER  ONE  OF  THE 
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LIT 

UB 

UM 

US 

FOALPH 

FO BETA 

FRALPH 

FRBETA 

HALPH 

HBETA 

DELTA 

ICYCLE 

AIC 

N 

THE  PROGRAM  WILL  LOOK  FOR  THE  VARIABLE  NAME  ENTERED  AND , 

IF  THE  VARIABLE  NAME  IS  NOT  FOUND,  THE  PROGRAM  WILL  RETURN 
ASKING  IF  ADDITIONAL  RE-INITIALIZATIONS  ARE  WANTED  (REFERBACK  8). 
IF  THAT  VARIABLE  NAME  IS  FOUND,  THE  PROGRAM  WILL  RESPOND  WITH: 

ENTER  THE  NEW  VALUE..: 

THE  USER  WILL  ENTER  THE  NEW  VALUE , AND  THAT  VALUE  WILL  BE 
ASSIGNED  TO  THE  VARIABLE  NAME  PREVIOUSLY  ENTERED. 

IF  THE  USER  ENTERS  ICYCLE,  THE  PROGRAM  WILL  RESPOND  WITH 
QUESTIONS  ABOUT  THE  RE-INITIALIZATION  OF  AIC  (REFERBACK  A) 

AFTER  CHANGING  THE  CURRENT  VALUE  OF  ICYCLE.  THE  NUMBER  OF 
ELEMENTS  IN  THE  ARRAY  AIC  DEPENDS  UPON  THE  VALUE  OF  ICYCLE. 

BY  RE-INITIALIZING  AIC  EACH  TIME  THE  VALUE  OF  ICYCLE  IS  CHANGED, 
THE  PROGRAM  INSURES  THAT  THE  CORRECT  NUMBER  OF  ELEMENTS  IN 
AIC  EXISTS.  IF  THE  USER  DOES  NOT  WISH  TO  CHANGE  A 
VARIABLE  NAME,  THE  USER  MAY  ENTER.4  NONE 

IN  ADDITION,  IF  THE  USER  FEELS  ‘CAUGHT’  IN  THE  RE- INITIALIZATION 
PROCESS  AND  WANTS  TO  RESUME  PROCESSING  OR  WANTS  TO  SEE  A LISTING 
OF  THE  INPUT  DATA,  THE  USER  SHOULD  ENTER:  • NONE 
AT  THIS  POINT.  THE  RE-INITIALIZATION  PROCESS  CAN  ALWAYS  BE 
RE-ESTABLISHED  BY  ENTERING J NO 

WHEN  THE  PROGRAM  ASKS  IF  THE  INPUT  DATA  IS  CORRECT  (REFERBACK  5). 


THE  USER  WILL  BE  ABLE  TO  CHANGE  THE  VALUE  OF  ONLY  ONE  VARIABLE 
(BY  ENTERING  VARIABLE  NAME  AND  VALUE  IN  TWO  ENTRIES)  AND  THE 
PROGRAM  WILL  ASK: 

(REFERBACK  8>: 

DO  YOU  WANT  TO  RE-INITIAL JZE  ADDITIONAL  VALUES..: 


IF 


.HE  USER 
PROGRAM 


i»  . 


WILL  RESPOND  AS  IN  RE- INITIAL IZATIOK 


PAG£  15  4UAL11X  rtmiXlUABIM 
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(REFERBACK  e>.  IF  THE  USER  ENTERS  NO  OR  Nr 

THE  PROGRAM  WILL  RESPOND  UITH  A LISTING  OF  INPUT  DATA 

FOR  USER  APPROVAL  <REFERBACK  5). 

(NOTE:  AFTER  TERMINATING  PROCESSING  OF  PCRC»  IF  THE  USER 

IS  PROCESSING  ADDITIONAL  PROGRAMS  BEFORE  LOGOUT r THE  USER 
MUST  ENTER.*  CONNECT  OUTPUT 

AFTER  A COMMAND  REQUEST ♦ OUTPUT  WAS  RETURNED  (DISCONNECTED) 
INITIALLY  TO  BYPASS  PRINTING  OF  NUMEROUS  TRIVIAL  ERROR 
WARNINGS. ) 
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COMMAND-  ATTACH fFILE2,PCRC 
PF  CYCLE  NO.  = 001 

COMMAND-  ATTACHfLIBfIMSLf ID=LIBRARYfSN=ASD 

FF  CYCLE  NO.  = 99? 

COMMAND-  LIBRARYfLIB 

COMMAND-  RETURN  OUTPUT 

COMMAND-  CONNECT  TAPE? 

COMMAND-  FTNfI=FILE2fL=0 


2.173  CP  SECONDS  COMPILATION  TIME 
COMMAND-  CONNECT  TAPE1 


COMMAND-  LGO 

ENTER  UTfUBfUMfUS  . . J 
ENTER  FOALPHfFOBETA  . .J 
ENTER  FRALPHfFRBETA  . .J 
ENTER  HALPHfHBETA  . .J 
ENTER  DELTA fICYCLEfN  ...* 

ENTER  AA ( I ) f 1 = 1 f 4 ..J  0 f . 016 f . 032 f . 064 

ENTER  COST ( J ) fJ=1 f4  ..J  3f3f3f100 

DO  YOU  WANT  MULTIPLE  VALUES  FOR  AIC..J  N 
ENTER  VALUE  FOR  AIC..J  5000 


2000f .2614E-03f1.00f3.00 
2f .0190 
2f .0190 

If  .1 

. 004  f 6 f 1900 
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CRACK  SIZE  DATA 


INITIAL  FLAW  SIZE  DISTRIBUTION— UEIBULL 


SHAPE  PARAMETER 

= 

FOALPH 

= 

2.0' 

SCALE  PARAMETER- 

= 

FOBETA 

= 

.0190 

CRACK  GROWTH  PARAMETERS 
INSPECTION  INTERVAL 

UT 

2000. 

GROWTH  RATE 

= 

UB 

= 

♦ 2614E-03 

USAGE  BIAS 

25 

UM 

= 

1.00 

FLT  BY  FLT  VARIABILITY 

US 

5= 

3.00 

PROBABILITY  OF  CRACK  DETECTION— WEIBULL 

SHAPE  PARAMETER  = HALPH  = 1.0 

SCALE  PARAMETER  = HBETA  » .1000 

REPAIRED  FLAW  SIZE  DISTRIBUTION— WEIBULL 
SHAPE  PARAMETER  = FRALPH  = 2.0 

SCALE  PARAMETER  = FRBETA  = .0190 

PROGRAM  PARAMETERS 

CRACK  LENGTH  INCREMENTS  = DELTA  * .004 

NUMBER  OF  INSPECTIONS  = ICYCLE  = 6 

INSPECTION  AND  REPAIR  PARAMETERS 

NUMBER  OF  ITEMS  = N = 1000 

INSPECTION  FLAW  SIZE  = AA(I)»I=1»4  = O.OOOr  .OlAr  .032»  .064r 
REPAIR  COSTS  = COST< J> t J=1 r4  = 3.00»  3.00r  • 3,00.  100.0 

0» 

INSPECTION  COST  = AIC ( I I ) t I 1=1 t ICYCLE  = 5000. t 5000. r 5000. » 

5000. t 5000. t 5000. » 


IS  CRACK  SIZE  DATA  CORRECT..:  Y 


0 


A 


♦ 

.4000000E-02 
.8000000E-02 
. 1200000E-01 
. 1600000E-01 
.2000000E-01 
« 2400000E-01 
« 2G00000E-01 
. 3200000E-01 
.3600000E-01 
.4000000E-01 
. 4400000E-01 
« 4800000E-01 
< 5200000E-01 
. 5600000E-01 
* 6000000E-01 
. 6400000E-01 
. 6800000E-01 
. 7200000E-01 
. 7600000E-01 
. 8000000E-01 
• 8400000E-01 
.8800000E-01 
• 9200000E-01 
« 9600000E-01 
. 1000000E+00 
. 1040000E+00 
. 1080000E+00 
. 1120000E+00 
. 1160000E+00 
. 1200000E+00 
. 1240000E+00 
. 1280000E+00 
MIN.  A 0.000 


F-  ( A ) 

0. 

. 7665424E+01 
. 1462606E+02 
. 2028449E+02 
• 2423527E+02 
. 2A31019E+02 
. 2657786E+02 
. 2530251E+02 

• 22875&7E  f 02 
. 1973837E  >02 
. 1631065E+02 
. 1294020E+02 
• 9875126E+01 
. 7259753E+01 
. 5147545E+01 
. 3523800E+01 
. 2330905E+01 
. 1490947E+01 
. 9228043E+00 
. 5530006E+00 
. 3210313E+00 
. 1806324E+00 
.9855436E-01 

• 521661  IE-01 
♦2679917E-01 
. 13367f/8E-01 
. 6476784E-02 
. 3049373E-02 
. 1395542E-02 

• 621 1 453E-03 
. 2689419E-03 
. 1 133165E-03 
. 4642349E-04 

MAX.  A .128 


F+<  A) 

0. 

. 1242439E+02 
. 2236075E+02 
» 286381 1E+02 
. 3106173E+02 
. 3027301E+02 
. 273421 1E+02 
. 2334291E+02 
. 1909167E+02 
. 1508462E+02 
. 1 15&215E+02 
. 8606656E+01 
. 621 7902E+01 
. 4354470E+01 
• 2952821E+01 
. 1937605E+01 
. 1230068E+01 
. 755b850E+00 
. 4492323E+00 
. 258631 1E+00 
. 1442508E+00 
• 7798125E-01 
. 4087872E-01 
• 2078919E-01 
. 1026121E-01 
. 491 7659E-02 
. 2289250E-02 
. 1035553E-02 
• 4553371E-03 
. 1947205E-03 
. 8100373E-04 
. 3279202E-04 
. 1290746E-04 


J 

LOW  - UP 

P(If J) 

1 

0.000-  .016 

.021678 

2 

.016-  .032 

.086637 

3 

.032-  .064 

. 120660 

4 

.064-  .128 

.009682 

K<  1)  = .238658 

E(C)  S(C) 

6655.  313. 


CUMULATIVE  TOTALS 
E ( C ) S(C) 

6655.  313. 


J<3 
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PROGRAM  LISTING 


too-  PROGRAM  PCRC< INPUT, OUTPUT, TAPEl-lOO, TAPE9-100) 

110=  COMMON/ 1NT/  A,NN,AK,XX<503) ,FMI <303) * FPL (503) 

120*  COMMON/COM/  CC<500. 3) ,NIP 

1 30=  COMMON/UCOM/  UT.UB.UM.UB 

140-  COMMON/FOCOM/  FOALPH. FOBETA 

130-  COMMON/HCOM/  HALPH.HPETA 

140-  COMMON/FRCOM/  FRALPHt FRBETA 

170-  DIMENSION  BPAR< 4 > ,AA< 5 ) , HH( 503 ) .PQ< 4 ) .COST ( 4) , AIC < 20 ) , XYZ < 3 ) 

100-  EXTERNAL  FMINU9 

190-  DATA  YES/ 3H YES/ , Y/1HY/ 

200-  DATA  BPAR/4*0.0/ 

210-1001  WC1TE  <9*1000) 

220-1000  FORMAT < 4X , *ENTER  UT*UB»UM*US  ..1*,4X> 

230-  READt.UT  *UB*UM*U9 

240-  WRITE  ( 9,1010 ) 

250-1010  FORMAT <4X,*ENTER  FOALPH , FOBET A . .t*,4X> 

240=  READ*  * FOALPH* FOBET A 

270-  WRITE  <9.1020) 

280=1020  F0RMAT<4X.*FNTER  FRALPH . FRBETA  . .t*,4X> 

290-  READ*. FRALPH. FRBETA 

300-  WRITE  <9.1030) 

310=1030  rORMATf 4X.*ENTER  MALPH.HBETA  ,»1*.4X> 

320=  RFAD*. MALPH.HBETA 

330-  WRITE  <9.1040) 

340=1040  F ORMAT  < 4X » *ENTER  DELTA. ICYCLE.N  . .t*.4X) 

350-  READ*. DELTA. ICYCLE.N 

340=  WRITE  <9,1030) 

370=1050  FORMAT  ( 4X  » *F.NTER  AA<I).I-1,4  ..I*,4X) 

380=  READ*. <AA<I), 1-1,4) 

390=  WRITE  <9,1040) 

400” 1 040  M1RMA T < 4X • 9ENTER  COS T < J > , J- 1 , 4 ..t*,4X> 

410=  READ* , < COB  T<J).J-1,4) 

420=  WRITE  <9, 0743) 

430-0743  FORMAT  <4X.*D0  YOU  WANT  MULTIPLE  VALUFS  FOR  A!C..!*,4X) 

440=  RFAD  9001, REPLY 

450”  IT  (RFPLY  .EO.  YES  .OR.  REPLY  .EO.  Y ) GO  TO  233 

440-  WRITE  <9,7454) 

470=  READ* , AICVLU 

480”  DO  101  I 1-1 , ICYCLE 

490=  AIC  f 1 1 ) -AICVLU 

300=101  CONTINUE 

510=  GO  TO  1500 

520-235  WRITE  <9,4543) 

530-4543  FORMAT  <4X.*ENT£R  VALUES  FOR  AIC--ONE  VALUE  PER  ICYCLE. . I*. 4X> 

540-  READ*. <AIC< II ) ,11-1, ICYCLE) 

550=1500  WRITE  <9,t) 

540=1  FORMAT  </////, 23X , *CRACK  SIZE  DATA*,/) 

570-2  FORMAT  <4X,*INITIAL  FLAW  SIZE  DISTRIBUTION— WEIBULL*. ) 

580-3  FORMAT  <0X.*SHAPE  PARAMETER  - FOALPH  « *,P4.1) 

590=4  FORMAT  <0X,*SCALE  PARAMETER  * FOBETA  « *,F4.4) 

400=5  FORMAT  </.4X,*CRACK  GROWTH  PARAMETERS*. ) 

610=6  FORMAT <0X,* INSPECT  ION  INTERVAL  * UT  - *,F7.0) 

620” 7 FORMAT  <0X,*GROWTM  RATE  - UP  " *,F10.4) 

430=0  FORMA f <0X,*IISAQE  BIAS  - UM  - *.F3.2) 

640-9  1 ORMAT  <0X.*FLT  BY  FLT  VARIABILITY  - US  - *,F3.2) 

450-  WF<  I f E < 9 » ? ) 

640=10  H1RMA1  < / . 4X  , *PROBAML  I TY  OF  CRACK  DF  TEC  I ION  -WE1BIILI  * , ) 

470  11  TORMAI  (0X,*SHAPE  PARAMETER  - MAI PM  t.F 4.1) 

400  1?  IOPMAI  <HX,*nCAI.F  F'ARAME  IFR  - HBF1A  = *,rA.4) 

490-  WR1  IF  <9f.«)F0AI  PM 

ZOOM  TORMAI  </.4X,*RFPAIRfD  FLAW  SI7F  D 1 0 T R I BU 1 1 1 IN  - WF  I Bl  II  I.  * . ) 
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710  = 14  FORMAT  < 8 X,* SHAPE  PARAMETER  = FRALPH  = *,F4.1> 

720-15  FORMAT  <8X,*SCALE  PARAMETER  * FRBETA  = *,F 6.4) 

730=16  FORMAT  ( /, 4X » *INSPECT ION  AND  REPAIR  PARAMETERS* r ) 

740= 1 7 FORMAT  < 8X * *INSPECT ION  FLAW  SIZE  - AA<I), 1-1,4  = *, 

750=  • 4<F5.3,*,*)>  * 

760=18  FORMAT  <8X,*REPAIR  COSTS  * COST< J) , J-l » 4 * *f 
770*  ♦ 4<F8.2»*,*>> 

780=  WRITE  <9,4)r0BETA 

790” 19  FORMAT  < 8X , * INSPECT  ION  COST  - AIC<II),II=1» I CYCLE  * *, 

800-  • 10<r7.0,*,*)/8X»2<F7.0,*,t) ) 

010-20  FORMAT  < 8X  * tNUMDER  OF  ITEMS  - N = *.I6) 

020*21  FORMAT  < / , 4X , 4FR0GRAM  PARAMETERS** ) 

030=22  FORMAT  (8X»*CRACK  LENGTH  INCREMENTS  - DELTA  * *,F5.3) 

040*23  FORMAT  < BX * *NUMDER  OF  INSPECTIONS  - ICYCLE  * *»I2,> 

850=24  FORMAT  </////, 4X , *IS  CRACK  SIZE  DATA  CORRECT • * 1 $ » 4X ) 

860*  WRITE  <?»«> 

870*  WRITE  < 9 • 6 )UT 

fino-  WRf  IF.  <9*  7 )OB 

090  * URfTE  < 9 * 8 ) I IM 

700  WRITE  <9*9  MIS 

910*  WRITE  <9,10) 

920  WRITE  <9,11 )HALPH 

930*  WRITE  <9* 12)MRETA 

940-  WRITE  <9,13) 

950=  WRITE  <9,14 )FRALPH 

960*  WRITE  < 9 * 15 1FRBETA 

970*  WRITE  <9,21) 

980*  WRITE  <9,22)DELTA 

990*  WRITE  (9,23) ICYCLE 

1000=  WRITE  <9, 16) 

1010*  WRITE  (9,20)N 

1020=  WRITE  (9,17) AA < 1 ) * AA ( 2 ) , AA ( 3 ) , AA < 4 ) 

1030*  WRITE  <9, 18  >COST ( 1 ) ,COST (2) ,COST (3) ,COST  < 4 ) 

1040*  WRITE  <9, 19)<AIC<II), 11*1, ICYCLE) 

1050=  WRITE  (9,24) 

1060*  READ  9081, REPLY 

1070*  IF  < RF.PL  Y .NE.  YES  .AND.  REPLY  .NE.  Y > GO  TO  6001 

1000”  ECT0T*0.0 

1090=  SSOTOT-O.O 

MOO-  DO  100  1 1*1,  ICYCLE 

1110-  NN-II 

I 120-  IT  < II  .01. 1 ) CALL  IC8ICU<XX,FrL,NIP,WAR,CC»NIP-l , 1ER  ) 

II  10*  DO  40  1*1 ,501 

1 I 40*  A-TiFI  IA*FLOAT<  I 1 ) 

IIVl*  XX  < I ) -A 

1160-  Ml* I 

M 70  AA<5)*-A 

I I no  * FMT  < I >~UCAfiRF<FMIMUS,  -3. 0,3.0,  1 .F  4, 1 .T  A ,F  RROR,  HR) 

1190*  MM<  M*M<A)*FMI<I) 

1200-  If  <!  .F(i.  1)  GO  TO  40 

1210*  IlMIMKI)  rMMl-m.LT.  0.0  .AND.  EMI  < I ) .1.1  . l.OE-4)  GO  10  45 

1220*40  CON  T INI  If 

1230*45  NTT-N1 

1240  AK-0.0 

1250”  CAM.  ICSICU<XX,HH, NIP, SPAR, CC, NIP-1, IFR) 

1260*  DO  50  1*1,4 

1270-  DO  1313  IXYZMvS 

1200  - XYZ<  tXYZ  )*- AA<  IXYZ  ) 

1290-  IT  < XY7< IXYZ)  .GT . AA<5) ) XYZ(IXYZ)  = AA<5) 

1300**1313  CONTINUE 
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1 310^  CALL  PCSODIM  XX.HM.NIP.CC.NIP-1  »XYZ<  I )»XYZ(I+1)  »PQ<  I ) . IER) 

I 320-  AK=AK*PO<I> 

1 J30-50  CONTINUE 
1340=*  DO  60  I“t  *HIP 

1 330-  A = XX  < I > 

1360=  FPL  ( I )«FFLUS<A»  I ) 

1370-60  CONTINUE 
1 500~  WRITE  <9.2000)11 

1 190=2000  F ORMA  T < INI » ///4X » 20  < 1H* ) » 1 2 » * TH  CYCLE  *.20<1H*)) 

1400  WRITE <1.2005) 

1410=2005  F0RMAT<///8X.1HA.12X.*F-<A)*.12X.*FF<A)*./) 

1420=  HO  65  1=1. NIP 

1410=  UR I TE< 1.2010)  XX < I ) .FMI < I ) .FPL < I > 

1 440=2010  FORMAT < 1 X . 3 < E15 . 7 . IX ) ) 

1450=65  CONTINUE 

1460-  WRITE  <9. 201 2 ) AA<  5 ) 

1 470=201 2 r ORMAT <4X»*MIN.  A 0.000  MAX.  A *.F5.3> 

1480-  WRITE  <9,2015) 

1490  7015  FORMA  T <///l IX. 1HJ » * LOW  - UP* . 6X » *P< I . J ) * ) 

1500=  HO  70  1=1,4 

1510  WRITE  <9.2020)  I ,XYZ< I ) .XYZ< 141 ) ,FQ< I ) 

1520*2020  r ORMAT < 1 IX. II , 2X ,F5 . 3 , *-* ,F5 . 3 , 2X ,F1? , 6 ) 

1530=70  CONTINUE 

1540-  UR IFF  <9.?030)  II.AK 

1550*2030  F ORMAT <//12X.*K<*. I2»* )*. 3X. *»*. 3X.F12.6 > 

1560  - EI>A1C  < 1 1 ) 

1570=  SSOC’O.O 

1500--  no  00  j-l.4 

1590-  EC  *EC I FLOAT  <N)*PO<  J)*COST  < J) 

1600  - SSnC=SSOCfCOST< J)**2*FL0AT<N)*PQ< J)*<1 .0-P0<J) ) 

1610 "80  CONTINUE 

1620=  WRITE  <9.2040)EC»  SORT  < SSQC ) 

1630  -2040  T ORMAT <//13X.*E<C>*»9X»4H9<C>  . / . 6X . 2<  3X . F 1 2 . 0 ) ) 

1640=  EL  1 Of  *FC  TOT  TEC 

1650=  SSOTOT-SSOTOTTSSaC 

1660*  WRITE  < 9.2045)ECT0T.S0RT<SS0T0T) 

1670=2045  rURMAT  <//»  1 4X  . *CUMl)L  AT  I VE  TOTALS*  ./ 1 3X  . *E  < C )*.  9X . 4MS  < C ) . 

1600=  t /6X»2<3X»E12.0>> 

1 690*100  CONTINUE 

1 700=C  WRITE  < 9 » 2050 ) EC TOT  * SORT  < SSOTOT ) 

1710  *7050  FORMAT  <IM1 »RX»10<1H*)»*  TOTALS  *. 10 < 1H* ).///. 10X . *E <C ) -TOTAL* . 
1720-  « 4X . 10HS  <C)“T0TAL»/»6X»2<3X»E12.0) ) 

1730  9881  r ORMAT  < A10 ) 

1 740" 9876  F ORMA T<4X»*ENTER  THE  NEW  VALUE . . S*.4X) 

1 750 1 WRITE  <9.6000) 

1760-6000  rOPMAT  </////, 4X . *P0  YOU  WANT  TO  CONTINUE  PROCESSING. • J*»4X> 

1770=  REAP  9881  .RFPl.Y 

1700'  IT  < REPLY  .NE.  YES  .AND.  REPLY  .NE.  Y ) GO  TO  200 

1/90-6001  WRITE  <9.7000) 

IGOO  7000  I (IRMA  T <4X,*P0  YOU  UANI  TO  Rf-lNlUAtl/E  All  VAl  UF8 . . 1 * . 4X  > 

1810=  RFAP  98R1.RFP1Y 

1820-  If  < RF  PL  T .TO.  YES  .OR.  REPLY  .FO.  Y ) 00  TO  1001 

11130  UR!  If  <V . 77  77) 

1840'  REAP  Vfini.RIFlY 

1850=  rr  <Rf  fi  Y .NT.  YES  .AND.  RE  PL  Y . Nl  . Y > GO  TO  7002 

I860-  WRTTT  <9.1050) 

1070=  REAP*. <AA<  f ) . 1 = 1 .4  ) 

1800  7007  UR  I IF  (9.  7776) 

1TIV0-  7776  I ORMAT  <4X.*P()  YOU  WANT  TO  RE  TNJTIAtl/F  IMF  VAl  UE  OF  COS  T < J ) . . I * . 
1900-  I4X ) 
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1V10 

kl  AI»  V8H 1 » KEPI  Y 

1920 

11  (KEPI  Y . NL  • 

YES  . AND 

. REPLY  .NE.  Y ) GO  TO  7001 

IV  JO 

UK  4 It  < V » 1 060 ) 

1 9*0 

read*.  <i:osr<  J)» 

J I » 4 ) 

1950*  '001 

wkirt  tvroooo) 

I960  UOOO 

nilcM.11  C.X, .ENTER  IMF  VARIAUIt  NAME  UHUb'E  VALUE  IS  1(1 

1 V /<> 

•i  mu. . :*»4x> 

1 VHO  - 

MAH  988 1,  RE  PI  Y 

1 WO 

IF  (MEM  Y .Nl  . 

2HUT ) GO 

TO  8001 

2000  • 

UK  II I < V • 9876 ) 

2010- 

Kt  AH*. Ill 

2020- 

GO  IU  9VVVV 

20 JO- HOO 1 

U < KEPI  Y .NE. 

2MIIB)  GO 

ro  V001 

2040 

UK  lit  <9,VU76> 

7050- 

READ*,llb 

2060“ 

<»0  10  9VV99 

20/0  //// 

FUkhAI  <4X,*DU 

YOU  UANT  TO  RE- INI T I AL I ZE  THE 

VALUE 

2080 

• » 4X  ) 

20 VO “VOOl 

IF  < REPLY  .NE. 

2MUM)  00  TO  10001 

2100- 

UR1IE  < 9,9876 ) 

21  10- 

READ* ,UH 

2120- 

G(J  TO  99999 

2130=10001 

IF  (REPLY  .NE. 

2HUS)  GO  TO  11001 

2140 

MkITE  <9»90/6> 

2150= 

RE AH*. US 

2 1 60= 

GO  TO  99999 

21  70- 1 1 001 

IF  (REPLY  .NE. 

6MF0ALPH)  GO  TO  12001 

2100  - 

URITE  (9,9876) 

2190  = 

READ* ,FOAl  PH 

2200= 

GO  TO  99999 

2210- 12001 

IF  (REPLY  .NE. 

6HFOI1ETA)  00  TO  13001 

2220“ 

URITE  (9,9076) 

7230  - 

READ*.FODFTA 

7240 

GO  TO  99999 

2250-1 JO01 

IF  (REPLY  .NE. 

6MF RALPH)  00  TO  14001 

2260 

URITE  <9.9876) 

27/0  * 

READ*,r RALPH 

2280- 

GO  10  99999 

7290  14001 

IF  (REPLY  .NE. 

6HFRBETA)  00  TO  15001 

7300 

UR  HE  (9,98/6) 

2310  = 

KL  AH* , FRRE TA 

7370 

GO  TO  V9999 

7330  15001 

IF  (REP1Y  .NE. 

5HMAIPH)  00  TO  16001 

23  40^. 

URITE  (9,98/6) 

2 450 - 

READ*, MAI  PM 

2360- 

GO  TO  99999 

73/0  -16001 

IF  (REPLY  .NE. 

5HHDETA)  QO  TO  17001 

7300  - 

URITE  (9,9876) 

7390 

READ* , HUE  T A 

2400 

GO  1U  99999 

24 J 0-1/001 

IF  (REPLY  .NE. 

5HDELTA)  00  TO  18001 

2 420- 

URITE  (9,9876) 

2430  - 

READ*. DELTA 

74  40=* 

GO  TO  9999V 

2450=10001 

IF  (REPLY  .NE. 

6HICYCLE ) 00  TO  19001 

2460  j 

URITE  (9,9876) 

2470* 

Rl  AD* » ICYCl.E 

2400“ 

UR HI  <9,8765 

) 

2490- 

READ  9BB1. REPLY 

2500- 

If  (REPIY  .NE. 

YES  .AND.  REPLY  .NE.  Y ) GO  TO 

8642 
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IMN)  GO  TO  1 8003 


L'MO* 

2520- 1V001 
7530* 

ITjTO3 

7550 

2560*18003 

25/0 

2500* 

2590 

2600-7654 
76 1 0 : 0642 


GO  TO  236 
If  (RCPIY  .NE . 

WRITE  <9,9876) 

KT  All*  , N 
00  TO  9999V 
1»  (REPIY  • Nt  • 3HAIC)  00  10  99999 
WRITE  (9,8765) 

RE  AH  980 1 , REPl  Y 

U <Rk.Fl  Y .EU.  YE5  .OR.  REPLY  .10.  Y ) GO  TO  236 
FORMAT  < 4X , *E  N TER  VALUE  FOR  AIC..J*,4X> 

WRITE  < 9,7654) 


?6?0  - 
76  50--- 
2640- 


2650- 207 


2670-236 


RE  All*  rAICVI  U 
HO  202  1 1*1 , ICYCl E 
All  < II  )=AICVLU 
CONTINUE 
GO  10  99999 
WRITE  (9,6543) 


7680-  RE  AIi«,  <AIC<  II  ) , 11*1  , ICYCLE) 

7690  99999  CONTINUE 

2700-  WRITE  (9,1234) 

2710-1234  f UKrtA r (4X..H0  YOU  WANT  TO  RE-INITIAL  WE  ADDITIONAL  VALUES. 
2720-  READ  V081. REPLY 

-7J°-  IE  (REPLY  .EQ.  YES  .OR.  REPLY  .Ell.  Y ) 00  TO  7001 

2740-  GO  TU  1 GOO 

2/S0-200  CONTINUE 

2 7 AO---  STOP 

2770-  ENP 


1».4X> 


2 /HO- 

2 no 
2000- 
2810- 
2020- 
2b  30 
2840- 
2050- 10 
2060-20 
2870- 
2080 
2890- 
2900  - 


2910’ 
2920  = 
2930* 
2940“ 
2950= 


2960- 

29 /0  - 

2980- 

2990- 

3000- 

5010= 

3020--*  10 

3030-20 

3040- 

3050* 


3060" 

3070“ 

3000= 

3090* 

3100" 

3110-10 

31 20-70 

3130" 

J140- 


f UNCI  Kir#  F N I NUS  < X ) 

CUNNUN/INI/  A.  1 1 » AR  t XX < 503 ) »F  hi  < 503 >»F PL  (503) 
UP7U*U2D<X> 

IIP2  A<MP2D 

IFUI  .1*1.  1)  UU  IU  10 
T ENP-F  0 < Ul'2  ) 

(ill  IU  20 
UNP-FPI  I N ( UP 2 ) 

IENF’2-G<  X ) 

I NINUS=UP2D<TENP<TENP2 

IF  < FHINUS  .LT.  0.0)FNINUS=0.0 

RETURN 

END 


FUNCTION  U2D( X ) 

CONNON/UCON/  UTfUBfUHfUS 

U2D"EXP ( -UB< ( UN<UT  4SQRT < UT ) <USS<X ) ) 

RETURN 

END 


FUNCIION  FO(X) 

(QNNON/FOCON/  F OALPH r FOBETA 
IF  <X  .LE.  0. ) UO  TO  10 
TERN1=<  X/FOBE  TA  > <<  < FOAl.PH-  1 • ) 
TERN2-  < X/FOBE TA ) <<F  OALPH 
GO  TO  20 
TERN1  =■  TERN2-0 » 

FO-FOALPH/FOBE TA<TERN1<EXP<  TERN2) 

RE  1 URN 

END 


FUNCTION  M<  X > 
CONMON/HCOH/  HALPH*HBETA 
IT  <X  .LE.  0.)  GO  TO  10 
TERN—  < X/HBETA  ) «HALPH 
00  TO  20 
TFRN-0 • 

M" 1 .0  FXr<TERN) 

Rr 1 URN 
END 
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mint:  r i on  fMua<x,  i> 

< <IH«0N/1NI/  A » i 7 . AK  y»/4  A', 

ht turn <<*K*i ***** m X > ; *r«i i'jf ,* 1 505 ' 

END 


3200= 
3210  = 
3220= 
3230  = 
3240- 
3250= 

32 AO- 10 
3270=20 
32U0- 
3290  = 


FUNCTION  FR<X> 

CONMON/FRCOM/  FRALRH,FRBE1A 

TERM 1 = <X/FRPFT  < pr^LRH- 1 ) 

GoRro*;oX/FRPETA>**FRALF,M 

'E KHI-1ERM2*0. 

lEZ^LrH/FR,,E'a$1Er>M,$EXf<J^> 

END 


HiNL  I ION  0<  X ) 

^;uR°hsort<2-°*3-"i5’2«‘',«-‘-x**^.o» 

END 


FUNCTION  FPLIN(X) 

COMMON/COM/  CC  f 500^3 w NIP 3 ’ ’ FHI < 503  > ' FPL  < 503 ) 
7 < 1^E^X ION  2<1)rS0l-<l> 

.IC9fVU(XX»FPLpNIP»cr»NiP  « 7 cm 
Pr*LIN-SOI  <1  ) nirf  cc *NIP— 1 »Zp SOL * 1 » IER ) 

PF TURN 
END 
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